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This article first reviews the basic physics of rotating stars and their evolution. We examine in 
particular the changes of the mechanical and thermal equilibrium of rotating stars. An important 
(predicted and observed) effect is that rotating stars are hotter at the poles and cooler at the 
equator. We briefly discuss the mass loss by stellar winds, which are influenced by the anisotropic 
temperature distribution. These anisotropics in the interior are also driving circulation currents, 
which transports the chemical elements and the angular momentum in stars. Internal differential 
rotation, if present, creates instabilities and mixing, in particular the shear mixing, the horizontal 
turbulence and their interactions. 

A major check of the model predictions concerns the changes of the surface abundances, which 
are modified by mass loss in the very massive stars and by rotational mixing in O- and B-type 
stars. We show that the observations confirm the existence of rotational mixing, with much larger 
effects at lower metallicities. 

We discuss the predictions of stellar models concerning the evolution of the surface velocities, the 
evolutionary tracks in the HR diagram and lifetimes, the populations of blue, red supergiants 
and Wolf-Rayet stars, and the progenitors of type Ibc supernovae. We show, that in many 
aspects, rotating models provide a much better fit than non-rotating ones. Using the same physical 
ingredients as those which fit the best the observations of stars at near solar metallicities, we 
explore the consequences of rotating models for the status of Be stars, the progenitors of Gamma 
Ray Bursts, the evolution of Pop III stars and of very metal poor stars, the early chemical evolution 
of galaxies, the origin of the C-cnhanccd Metal Poor stars (CEMP) and of the chemical anomalies 
in globular clusters. 

Rotation, together with mass loss are two key physical ingredients shaping the evolution of massive 
stars during the whole cosmic history. 
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I. INTRODUCTION 

It may be surprising that we emphasize the need of accounting for the effects of rotation in stellar evolution, since 
in practice all stars are rotating around their axis. However, it must be recalled that the so-called "standard theory" 
of stellar evolution generally ignores the effects of stellar rotation and treats the stars as non-rotating bodies, despite 
the fact that stars more massive than about 1.5 M Q rotate fast on the average (Fig. [I}. 

The recent progresses in astrophysical observations, particularly in high resolution spectroscopy and in asteroseis- 
mology, have shown many significant deviations from the standard models, for example the many large nitrogen 
enrichments resulting from mixing in massive stars (Sect. IVII[) . These observations show the need to also account 
for the various effects of rotation in stellar modeling. All model outputs are finally modified by the proper account 
of rotation: the stellar luminosities and radii, the lifetimes, the chemical abundances at the surface, the helio- and 
asteroseismology responses, the nature of supernova explosions, the amounts of nucleosynthetic products, the nature 
of the final remnants, . . . 

For practical purpose, it is often convenient to distinguish 4 main groups of rotational effects in stellar physics. 

1. - The equilibrium configuration of rotating stars: it results from the centrifugal force on the stellar equilibrium. 
The equipotcntials arc modified and in particular the shape of the stellar surface, which affects the surface T e g and 
gravity distributions. 

2. - The effects of rotation on mass loss or accretion: in fast rotating stars, the isotropy of mass loss (or accretion 
when present) is destroyed and anisotropies appear and are effectively observed. 

3. - The rotational mixing: the internal distortion induces circulation currents which transport the elements and 
angular momentum, while differential rotation may produce several instabilities which also contribute to the transport 
processes. 

4. - The interactions with magnetic field: the presence of an internal magnetic field may produce an internal coupling 
of rotation, leading to solid body rotation, while external fields produce some magnetic braking. A major uncertainty 
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concerns the existence of a dynamo in radiative regions with differential rotation. We examine some properties of 
such a dynamo. 

The above distinction is evidently a simplification since the various effects are related. For example, it is the 
modification of the internal equilibrium structure which drives the mixing. In turn, the mixing modifies the internal 
distribution of the elements and this also influences the equilibrium structure. 

In this review, we focus on the rotational effects in the Main-Sequence (MS) and post-MS phases, i.e. in the nuclear 
phases. The effects of rotation in star formation are also a most important chapter of astrophysics, but they would 
deserve another specific review. We consider here the case of single stars, the many effects of rotation in relation with 
tidal interactions in binaries are also beyond the scope of this review. Most of the effects discussed here in the case 
of single stars evidently have their counterparts in binaries, however their modeling is still in its infancy. 

Recent reviews on the observational and theoretical aspects of stellar rotation are given in IAU symposi um 215 
(|Maeder fc Eenensl , |2004) and the theoretical aspects have recently extensively been reviewed (|Maederl . l2009f) . 



II. THE MECHANICAL AND THERMAL EQUILIBRIUM OF ROTATING STARS 

The equilibrium and stability configurations of rotating stars have been reviewed long ago ( Lebovitzl . Il967fl . In 
practice, except for stars with little internal density contrast, like white dwarfs or neutron stars, the approximation 
of the Roche model is acceptable. It assumes that the effective gravity results from the matter centrally condensed, 
with in addition the effect of the centrifugal force. For all stellar masses, the rotational energy of the Roche model 
represents at most about 1% of the absolute value of the potential energy. 

The properties of rotating stars depend on the distribution of the angular velocity f2(r) inside the stars. The 
simplest case is that of solid body rotation, i.e. f2 = const., while more elaborate models include differential rotation. 



A. The mechanical equilibrium for uniform rotation 

Let us first consider the case of a constant angular velocity Q in the Roche model in hydrostatic equilibrium. The 
gradient of pressure P is given by 

- VP = -V$ + - il 2 V(r sinz?) 2 , (1) 
Q 2 

g is the local density, $ is the gravitational potential, which is unmodified by rotation in the Roche approximation and 
which gives the gravitational acceleration g = — V$ = — G ^ r - , r being the distance to the center. The components 
of g are (—g, 0, 0) and g = If £1 is constant or has a cylindrical symmetry, the centrifugal acceleration can also 
be derived from a potential, say V. In case of constant f2, one has 

~ W = Q 2 m and thus V = -]^l 2 w 2 . (2) 

zd = r sin$ is the distance to the rotation axis. The total potential ^ is ^ = $ + V and with the usual Poisson 
equation one has V 2 ^ = V 2 $ + V 2 V with V 2< f> = AnG g. In a uniformly (or cylindrically) rotating star, the 
centrifugal force derives from a potential and one has {V 2 V) ZU = ^ (— vj 2 Q 2 ^ = — 2 The Poisson equation 
becomes 



and the equation of hydrostatic equilibrium 



V 2 * = A-kGq-IW . (3) 



-VP = -V* = gos- (4) 



The effective gravity g c s results from both gravitation and centrifugal acceleration (take care of the sign of $ and , E r ). 
These expressions imply that the pressure is constant on an equipotential, i.e. P = P(*£). The equipotentials and 
isobars coincide and the star is said barotropic, otherwise it is said baroclinic. Eq.Q shows that the density is also a 
function of *f? only. Through the equation of state P = P(g,T), one also has T — T(^>). The quantities g, P, T are 
constant on the equipotentials $ = const. 

The stellar surface is an equipotential. The total potential at a level r and at colatitude "d (t? = at the pole) is 

*( r> 0) = -O^L -\n 2 r 2 sin 2 § . (5) 
r 2 
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For a star of total mass M, let us call the stellar radius at colatitude Since the centrifugal force is zero at 

the pole, the potential at the stellar pole is just G M/R p , where R p is the polar radius. This fixes the constant value 
of the equipotential at the stellar surface, which is given by 



GM 1 , , , G M 



(6) 



This equation also applies to moderate differential rotation, if 51 can be considered as constant on an equipotential 
surface. 

The shape of a Roche model is illustrated in Fig. [2] for different rotation velocities. The effective gravity resulting 
from the gravitational potential and from the centrifugal force is given by (j4]). If e r and e# are the unity vectors in 
the radial and latitudinal directions, one has 



9cS 



GM 



+ Q 2 sin 2 ■& 



e r + [fl 2 i?(z9)sintfcostf] e# 



(7) 



Care has to be given for many applications that the gravity vector is not parallel to the vector radius, with an angle 
e given by cose = — f/cff • r/(\g c s\ ■ \r\) between the two vectors. 



B. Differential rotation: the case of shellular rotation 



The most studied case of differential rotation is that of the so-called "shellular" rotation ( Zahnl . Il992j ). i.e. with a 
rotation law f2 constant on isobaric shells and depending to the first order on the distance r to the stellar center, 



(8) 



where P 2 is the second Legendre polynomial. This rotation law results from the strong horizontal turbulence in 
differentially rotating stars, which imposes a constancy of on isobars, while in the vertical direction, the turbulence 
is weaker due to the stable density stratification. The above writing of fl(r, $) implies some limitations, since isobars 
are not identical to spherical surfaces, in particular it only applies to lo w or moderate rotation y elocities. This case 
is non-conservative. The surfaces defined by = const, are isobars ( Mevnet fc Maederl [T997). but they are not 
equipotentials and the star is baroclinic. 

The equations of stellar strucure for a star in shellular rotat ion can be written i n a fo rm close to that of the non- 
rotating case in order to minimize the necessary modifications (Me vnet fc Maederl Il997l ). One associates a radius rp 



to an isobar, it is defined by Vp = (Air/3) rp , where Vp is the volume inside the isobar. For any quantity like g e g, 
which is not constant over an isobaric surface, a mean value is defined by 



< 5cff > = 7!- f 9cs da , 

J*=const 

where Sp is the total surface of the isobar and da is an element da = 
Lagrangian coordinates (where Mp is the mass internal to an isobar), the equa tions of hydrostatic equil ibrium, of 
mass conservation, of energy production and of energy transport can be written (jMevnet fc Maederl . [l997t) . 



(9) 



sin dipd-d/ cos e of this surface. In 



dP 
dM P 



G M P 

"1 — r jp 

iirrn 



(10) 



drp 1 
dM P ~ Airrp-Q 



(11) 



dLp 
dM P 



£nucl 



■■grav j 



(12) 



dhxT 
dM P 



G Mp , 
-j -Jpmm 

Anro 



St 

V a d, Vrad-^— 
JP. 



(13) 
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Anrj, 1 f Airrp 



wMth / P = — — and f T = — ^ — . (14) 

GM P Sp < 9eff x > V 5 p / < 5off >< ff cff > 

where e„ MC / , e„ and e ffra D are respectively the rate of nuclear energy production, of neutrino losses and of gravitational 
energy production or absorption by mass unit. V a d and V ra( j are the adiabatic and radiative temperature T gradients 
(rflnT/dlnP). ~g is formally here the average density between two isobars. The above set of equations allows one to 
construct models of stars in differential rotation (shellular case) with little modifications with respect to the usual set 
of standard equations. The results are obtained as a function of Mp the mass of an isobar of average coordinate rp. 
In practice, rp is the radius of the isobar at a colatitude -d such that P2(cos$) = 0, i.e. w 54 degrees. 



C. The critical velocities 

The critical or break-up velocity is reached when the centrifugal force becomes equal to the gravitational attraction 
at the equator, i.e. when 

^crit = d3 > (15) 

-""e.crit 



where i? e ,crit is the equatorial radius at break-up. With f2 cr ;t in the equation of the surface © at break-up, one gets 
for tf 
thus 



for the ratio of the equatorial to the polar radius at critical velocity, ^ c ' crit t — ^. The equatorial break-up velocity is 



2 o2 p2 GM 2GM 

«crit = ^crit #e,crit = ~~5 = o~p • ( 16 ) 

-^c,crit 3 -ftpjCrit 



Introducing a non-dimensional rotation parameter ui, defined as the ratio of the angular velocity to the angular 

i city at 
x = R/R p ,, 



velocity at break-up, u> — i.e. uj 2 = — qm"* > ^ 2 = Jf anc ^ ^ e ec l ua tion of the surface © becomes with 



1 i 2 2 • 2 q ^?p,crit f-tn\ 

- + —uj z x z sm <& = . 17 

x 27 Rp(w) 

If the polar radius would not change with lj, the second term would be equal to 1; Eq. (jTTJ) is thus an algebraic 
equation of the 3 rd degree, which gives the values of the radius as a function of U) and I?. For low rotation (ui < 0.5), 
a linear approximation of v in term of uj is valid 



8 GMoj 2 „ 8 GMuj 2 4 „ ., 

v 2 = 5 Rl w = -uj 2 v 2 - t . (18) 

27 Rl CTit ° 27 i? p , crit " 



Fig. [3] shows the critical velocities v CT it for stars of various masses and metallicities. The critical velocities grow with 
stellar masses, because the stellar radii increase only slowly with stellar masses. The critical velocities are very large for 
low metallicity stars, since their radii are much smaller as a result of their lower opacities. In the most massive stars, 
where radiation pressure nearly compensates gravity, the effects of radiation modify the expressions of the critical 
velocity making them lower. In this case, one speaks of the fir limit, since the stars rea ch their critical velocity, d ue 
to the combined effects of high radiation (the Eddington limit) and of centrifugal force ( Maeder fc Mevnetl . l2000al) . 



D. The breakdown of thermal equilibrium in rotating stars 

Rotation deeply affects the internal thermal equilibrium. As seen above, the equipotentials are differently spaced 
as a function of the colatitude i?, since the star is distorted. The equipotentials are closer in polar regions and more 
separated near the equator. The radiative flux at some latitude is proportional to the local gradient between the 
equipotentials (the effective gravity, Sect. III.EI) . If there is no local energy source or sink, the total energy through an 
equipotential is conserved. Thus, there is (in a simple picture of uniformly rotating stars) an excess of flux along the 
polar axis and a deficiency near the equatorial plane. This thermal imbalance generates global circulation motions 
in the meridian plane. These currents then contribute to the mixing of chemical elements and transport of angular 
momentum (Sect. IIV|) . 



6 



E. The von Zeipel Theorem and observations 

The von Zeipel theorem relates the radiative flux at some colatitude i? on the surface of a rotating star and the 
local effective gravity g c ff (51, $), function of the angular velocity J! and i?. The flux is F(Cl, $) = — \ VT(f2, t?) with 
X = A< ^ c K T g ■ We consider the case of a uniformly rotating star. The equipotentials and isobars coincide, they are also 
surface of constant T and g. Thus, one has with @ 

F(Q, 0) = -x ^ VP(fi, 0) = -qx ^ 5cff (fi, 0) • (19) 

The pressure gradient and effective gravity are parallel. The term (gxdT/dP) is constant on a given equipotential, 
thus the flux is proportional to the effective gravity on the equipotential. This coefficient of proportionality can easily 
be estimated and one obtains 

= ~ 47r g M . 9esM (20) 
/ O 2 \ 

with M* = Mil — , (21) 

V 2irG g M J 

where ~g M is the average internal density. The above relation applies to any equipotential, but in general the von 
Zeipel theorem is considered at the surface of a star of total mass M and luminosity L. Eq. (f2"U]) is the von Zeipel 
theorem: it says that the radiative flux at the surface of a rotating star is proportional to the local effective gravity 
at the considered colatitude. From g c s , one obtains the radiative flux locally and thus T e g 



Both g c g and T e fi vary over the surface of a rotating star and influence the emergent spectrum. The equatorial 
regions are fainter and cooler than the polar ones, which are brighter and hotter (the differences of T c g may reach 
a factor of 2). This effect is c alled the gravi ty-darkening. The von Zeipel theorem in a star with shellular rotation 
shows only minor differences ( Maederl . Il999f) with respect to (f2"U)). The differences, which essentially depend on the 



fl gradient close to the surface, may slightly increase the contrast between the pole and equator. 2-D models have 
tested the validity of the von Zeipel relation ( Lovekin et all 120061 ). Depending on the rotation laws, there are some 



small differences. 

Interferometric observat i ons w ith the VLTI of the Be star Achernar indic ate a ratio R c /R p « 1.5 

( Domiciano de Souza et all |2003() . A recent analysis of the data by ICarciofi et al. (2008) confirms that the ob- 



servations well agree with a rigid ly rotating star at uj = 0.992 in the Roche model, provided it is surrounded by 
a small disk. VLTI observations ( Domiciano de Souza et all 120051 ) of the fast rotating star Altair (M w 1.8 Mq, 



A7IV-V) confirm a gravity-darkening as pred icted by the von Zeipel theorem. This is also supported by further 
studies for Altair ( Peterson et all [2008, 200 61), which rotates at 90% of its breakup angular velocity. The effect of 



gravity-darkening receives further support ( Monnier et all l2007t ). but instead of an exponent 0.25 as in Eq. (|2"2"j) . 



an exponent j3 = 0.19 is favored. Also, these authors notice an equatorial darkening stronger than predicted, which 
might result from faster equato rial rotation, of differe nces due to convection or opacity effects, etc. A value of /? as 
small as /3 — 0.08 is also found ([van Belle et al. , 2006) in the case of Alderamin (a Cep, type A7IV-V) , which rotates 



at 83% of its breakup velocity. Thus, the von Zeipel relation is confirmed, but there are some possible deviations. 



III. MASS LOSS FROM MASSIVE STARS 

Massive stars above 10 M Q lose, at solar metallicity, a significant amount of mass during their evolution. As an 
example, a star with an initial mass of 60 M loses about the half of its mass during its MS evolution. Due to mass 
loss in further stages, it is left with only about 10 M Q at the time of supernova explosion. The mass loss rates M of 
OB stars reach about 10 -5 Mq yr _1 with wind velocities up to 3000 km s _1 . For these stars M is a billion times 
larger than for the Sun, for which M is about 10 -14 Mq yr _1 .The mass loss results from the stellar winds driven 
by the strong radiation pressure of luminous stars, which pushes the mass outside. The absorption of radiation by 
the spectral lines is the main effect transferring momentum from radiation to matter. In red supergiants, mass loss 
also occurs due to the absorption and diffusion of radiation by dust. On the whole, mass loss is a dominant effect 
influencing all the outputs of stellar evolution and nucleosynthesis. 
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The physics of stellar winds has been reviewed in several studies ( Kudritzki fc Puis) . l2000t lLamers fe Cassinellil . 



Owing to the many uncertainties, stellar model makers generally apply expressions of M based on observations, 
which are also not free from problems. A major one is due to the clumping effects of matter in the wind, This has 
already lead to a reduction of the M by a factor 2 or 3 over recent years and questions arise about a further reduction . 
For hot stars, there are several parametrizations of the M rates deduced from theoretical models (|Kudritzki fc Pulsl . 
120001 : IVink et all . l2000l l200ll ) for different M, L, T cS and metallicities Z. For O stars, M (taken positive) behaves 
globally like 

M ~ L 16 and M ~ Z 7 . (23) 



Let us note that lMokiem et al] ( 20061 ) find that for stars with logL/L^ sup erior to about 5.4, the wind strengths arc 



in excellent agreement with the theoretical predictions of IVink et al. I (|200lh 



Dust-enshrouded re d giants and supergiants have strong mass loss rates due to the high dust opacity. 
van Loon et all l2005f) suggest the following empirical formula: 



logM = -5.65 + 1.05 log ( — ^— I - 6.3 log f I . (24) 



io 4 L Q y °V 3500 K 

Wolf-Rayet (WR) stars are bare cores left over by mass loss from massive stars. Their M are very high, due mostly 
to their hig h L/M ratios. Models indicate that the M values scale with the actual WR masses like ( Nugis fc Lamersl . 
120001 . l2002f> . 

logM = A + B\og— , (25) 



M 







with av erage values A = —5.73, B = 0.88 ±0.14 for WR stars in general. All these parametrizations are continuously 
revised ( Grafener fc Hamannl . 120081 ) . 



Rotation affects the mass loss in two ways ( Maeder fc Mevnetl . [2000b): - 1. it makes the stellar winds anisotropic 



because the polar regions are hotter so that the radiation pressure is stronger, - 2. for given luminosity L and T D g the 
average mass loss rates are increased. Fig. |4] illustrates the distribution of the mass loss around a very massive stars. 
The higher T in polar regions drives a stronger mass flux than in the eq uatorial regi ons. This characteristic peanut 
shape has been observed in several nebulae around very massive stars (|Notal . Il999h . For stars with an equatorial 



T e ff lower than about 23000 K, the equatorial opacity increases enough to allow the radiative acceleration to push 
(despite the lower flux) the matter outward and to produce an equatorial ring. Some stars, the so-called B[e] stars, 
simultaneously show a peanut-shaped nebula and an equatorial ring. 

The anisotropic mass loss rates may be integrated over the stellar surface to give the overall M value. The result 
is that if we consider two stars with the same L and T e g, b ut one with rotation veloci ty v and the other one with 
v = 0, the ratio of their M-rates may reach very high values ( Maeder fc Mevnetl . l2000bl ). especially for stars close to 



the Eddington limit (i. c. for the brightest stars where the radiation pressure gradient becomes close to gravity). For 
hot stars far enough from the Eddington limit, the effects are limited, amounting to an increase of at most 70 %. 



IV. THE CIRCULATION AND MAIN INSTABILITIES IN ROTATING STARS 

Rotation drives internal circulation currents and several kinds of turbulent motions may arise from differential 
rotation. All these motions may contribute to the internal transport of matter and angular momentum. 



A. The transports of elements and angular momentum 

Many effects of element transport in stars can be considered as a diffusion process characterized by a diffusion 
coefficient D. The kinetic theory of gases indicates that in an isotropic turbulent medium with an average velocity v 
and a mean free path I the diffusion coefficient is D = (1/3) £v. Let us call JQ, the abundance in mass fraction of 
particles i. If there is a gradient of composition VXi, the diffusion velocity of particles i is vl = —(D/Xi) VXj. The 
diffusion equation of the chemical elements can be derived from the continuity equation, this gives at a Lagrangian 
mass coordinate M r , 



^K^D^). (26) 



(/ r 2 dr \ dr 
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The boundary conditions at the center and surface are — 



Or 



M r =0 



0, where M is the total mass. The 



M r =M 



diffusion coefficients may vary a lot and great care has to be given to their interpolation, if necessary ( Mevnet et all 
12004ft . 

In a differentially rotating star, the evolution of the angular velocity fl has to be followed at each level r (for shellular 
rotation), so that a full description of f2(r, t) is available. The values of f2(r, t) influence the mixing of elements and 
in turn the evolution of 0(r, t) also depends on the mixing processes and on the distribution of the elements. The 
derivation of the equation for the transport of angular momentum is not straight forward. In the case of shellular 
rotation, the equation in the Lagrangian form becomes (jMaederl . 120091 : IZahnl . Il992t ) 



4^) 



1 d 



1 d 



— -( gr *nU 2 (r)) + -- gDr* — 



5 r 2 dr 



Or 



4 on 



dr 



(27) 



There, f2 is the average value of f2 on an isobar, or in approximation © the value at /-^(cos - )?) = 0. U 2 is the radial 
component of the meridional circulation velocity (see next subsection). The second term on the right is a diffusion 
term, similar in its form to (|26[) . while the first term on the right is an advective term, i.e. the transport by a velocity 
current. We notice that Eq. (|2"o| does not contain such an advective term. It could contain a term of that kind, 
however in general it is considered that the combined effect of turbulence and circulation currents is equivalent to a 
diffusion for the element transport (see Eq. I34p . 

Great care has to be brought to the solution of this equation, since U 2 given in the next subsection contains terms 
up to the third derivative of O, so that this equation is of the fourth order. If there is no viscous momentum at 
the edges, one has two boundary conditions given by ^ = at the limits. The two other conditions are given by 
integrating the transport equations at the bottom and top limits r^ and r t , 



I 



lgr 4 nu 2 

5 



gr i nU 2 



d 
df 



d 






n 


~ dt 




r R 


n / 


r 







m r = rb . 



+ Mn in r = r t 



(28) 



The radial components U 2 of the velocity at the surface and center are zero. A4q represents the momentum of force 
applied at the stellar surface, typically by magnetic field in solar type stars or by tidal effects in binary systems. 



B. The meridional circulation 

The breakdown of the thermal stability in rotating stars (Sect. III.D|) drives some circulation currents. For long there 
was a severe physical problem: the solutions for meridional circulation were not conserving the angular momentum. 
Thus the reality of the circulation was often questioned and solutions without circulation were envisaged. Zahn has 
made a great step forward by showing that one must treat simultaneously the equation for energy conservation which 
expresses the thermal imb alance, the Poisson equation and the conservation of angular momentum, in order to have 
a self-consistent solution (Zahi], I199S 



The derivation of the circulation velocity for shellular rotation is rather lengthy ( Maeder fc Zahd . ll998llZahnl . ll992h . 
The starting point is the equation of energy conservation, 

P T ^[=V- (xVT) + g e-V-F h , (29) 

where S is the entropy per unit of mass and \ the thermal conductivity (Sect. III.E[) . The term e refers to the nuclear 
energy production rate only. — — Dh gCp VhP is the flux of thermal energy due to the horizontal turbulence with 
a diffusion coefficient Z?h (Sect. IIV.C|) . Cp is the specific heat at constant pressure P and Vh expresses the horizontal 
gradient. On an isobar, the variables are expanded around their average up to the second Legendre polynomials 
P 2 (cosz9), where <d is the colatitude, for example T (P, z9) = T (P) + f (P) P 2 (costf). 

The velocity of meridional circulation is the main quantity characterizing this effect, it is expanded into 2 components 

, , x ^ / \ dP 2 (cos i?) _ , , 

U = U 2 (r) P 2 (cos t9) e r + V 2 (r) '- e# , (30) 
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Ui (r) is the amplitude of the radial component of the meridional circulation velocity, V2 (r) is the amplitude of the 
horizontal component. The general expression of U%(r) in a stationary situation is 



U 2 (r) 



P 



JgCpT (v ad -V + |V M ) 



L(r) 
M*(r) 



{Ec. 



(31) 



V M = (d\n fi/dlnP), ip = (d In g/d In aj)p,t, <5 = — (5 Ing/d hi T)p jfl , En contains dozens of terms up to the third 
derivative of il and E^ up to the second derivative of the mean molecular weight. An important term, also present 
in the case of uniform rotation in an homogeneous medium, is 



Er, 



1 



n 2 



2ttGq 



nv 

GM r 



(32) 



This term is proportional to the ratio of the centrifugal to the gravitational force. The second term in the square 
bracket is the so-called Gratton-Opik term. In the deep interior, it is generally negligible while it may become 
dominant in the outer layers, where the local average density ~g is negligible. The sign of U2 is however not determined 
only by the sign of the Gratton-Opik term and the complete expression has to be considered. Generally the inner 
cell is turning upward along the polar axis, while the external cell turns the opposite way (Fig. [5]). This means that 
the inner cell is bringing angular momentum inward, while the outer cell brings it outward. However, the circulation 
patterns, which depends on the derivatives of f2 and /z, may be different according to the stellar model considered. 

The circulation velocities are typically of the order of a few cm s^ 1 in upper MS stars. The velocities are larger 
at the beginning of the Main Sequence (MS) phase when rotation is (supposedly) uniform, then in a few percents 
of the MS timescale the velocities become smaller reaching an equilibrium. At the end of the MS phase, they may 
again become large (negatively) due to the very high Gratton-Opik term, because of the very low density in the outer 
layers. The circulation timescale t C irc, also called the Eddington-Sweet timescale t^s, is of the order of the stellar 
radius divided by the velocity, 



£es 



R 



GM 2 GM 



^KH 



U 2 {R) LR Q 2 R 3 Q 2 R 3 /(GM) 



(33) 



where <kh is the Kelvin-Helmholtz timescale, i.e. the thermal timescale of the star (about 3 x 10 7 yr for the Sun). The 
circulation timescale is equal to about £kh divided by the ratio of the centrifugal force to the gravity fi 2 i? 3 /(GM), 
thus £es is some multiple of £kh- 

Meridional circulation is important, because it transports both angular momentum and chemical elements. The 
presence of horizontal turbulence in differentially rotating stars affects these two kinds of transport differently. The 
interaction of the transport of chemical e lements by meridional circulation and horizontal turbulence results in a 
diffusion process (|Chabover fc Zahnl . Il992h . with a diffusion coefficient 



D fi 



30 D h 



(34) 



where -Dh is the diffusion coefficient by the horizontal turbulence as discussed in next Section. The above expression 
is valid only if -Dh is large with respect to r U2 and to the coefficient of vertical diffusion. 

The above reasoning does not apply to the transport of angular momentum, thus Eq. (|27p contains both advective 
and diffusive terms. The reason is that angular momentum conservation limits the effects of horizontal motions on 
the transport of the angular momentum. Indeed, turbulence may be present without a ne t flux of angular mom entum, 
however it maintains a shellular distribution of velocity, as shown by numerical models ([Charbonneaul . 119921 ). 



C. The horizontal turbulence in differentially rotating stars 

Rotation and especially differential rotation generates turbulent motions. On th e Eart h, we have the example of 
west winds and jet streams. In a radiative zone, the turbulence is stronger (IZahnl . 11992ft in the horizontal than in 



the vertical direction, because in the vertical direction the stable thermal gradient opposes a strong force to the fluid 
motions. The horizontal turbulence is characterized by a diffusion coefficient -Dh, which also exp resses the hori zontal 
viscosity. There are several observational indications supporting the existen ce of this turbulence (iMaedeit 20091). The 
expression of -Dh is uncertain and three different forms have been proposed ( Maederl . 120031 : iMathis et al. [ 12004 IZahnl . 



1992). That by Mathis et al. is based on the Couette-Taylor experiment, while the other two are phenomenological. 
The expressions by Mathis et al. and Maeder give similar results, larger by orders of magnitude with respect to those 
obtained with Zahn's expression. 
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A diffusion coefficient behaves iike £ 2 /tdis, where £ is a characteristic length (here typically r) and idiff a diffusion 
timescale. Thus, the question is to know which timescale is to be considered. One may see that the coefficient by 
Zahn corresponds to about a value tdis ~ r/Vz, where V2 is defined in Eq. (|30]) . where one has accounted that 
y 2 U2/3. The coefficient by Mathis et al. corresponds to a value t^iff ~ 1/(^2), where is defined in Eq. © 
and j3 is a numerical coefficient of about 1.5 • 10 -5 , determined empirically. The coefficient by Maeder corresponds to 
tdiff ~ [r/(Oa V2)] 1 ' 2 , meaning that the timescale considered is that necessary for the differential effects of the Coriolis 
force to cross the radius. At present, it is still uncertain which coefficient is to be chosen, even if the last two are in 
numerical agreement. 



D. Shear mixing and other instabilities 

A difference of horizontal velocity V between two layers may generate a shear instability. The vertical density 
stratification (along axis z or r) favors stability and the balance of the two effects has to be considered. The criterion 
expressing the instability is known as the Richardson criterion, it is 

_ g dg/dz 

Ki = - - - — 3 < Utah ■ (35) 

Q (dV/dzf 

The critical Richardson number is often taken equal to 1/4, some studies suggest a value up to 1.0 (jCanutol . l2002h . 
The criterion says that an instability develops if the excess of energy TZi cr i t g (SV) 2 in the differential motions (with 
respect to an average velocity) is bigger than the work g dg Sz necessary to exchange the matter vertically. The 
criterion can be expressed in term of the Brunt— Vaisala frequency N, i.e. the frequency of oscillatory motions under 
the restoring force of gravity, 

m = - N -g < fticrft , with N 2 = ^- ( V ad - V + ^V p ) . (36) 
{dV/dzf Hp \ 5 'V 

The density stratification, in particular in regions with a high /j, gradient, is generally able to prevent the shear 
instability. However, the thermal diffusivity of a displaced fluid element weakens the stabilizing effect of the thermal 
stratification and this may lead to a significant mixing. The left hand of Eq. (|36|) is then multiplied by v£/(6K), 
where v and I are the velocity and lengthscale of the vertical motions (the f actor 6 applies to a spherica l fluid element) 
K is the thermal diffusivity K = AacT^ / {2>ng 2 C p) . The criterion is 



(tne tactor o a pplies to a spnerica J 
then l|Maederl 11995 IZahnf . fl992l) 



N 2 vl 

K ' s m^™ <K "- (37) 

For constant /i, the coefficient for the vertical diffusion by shears is easily expressed. One supposes that the turbulent 
diffusion is dominated by the largest eddies satisfying the Richardson's criterion 

Ahear = ~V t = 2 TZi clit K } . (38) 

3 iv z 

This coeffcient expresses the transport of the elements and of the angular momen tum produced by shear turbulence 
in an homogeneous medium, with appropriate averages on isobars (jMaederl . 120031 . In principle, the thermal gradient 
is slightly modified by the turbulence due to the shears, however this effect is often ignored. 

In a zone with variable /i, the situation is different. The horizontal turbulence tends to reduce the composition 
difference between the turbulent eddies and their surroundings. T hus, the restoring for ce of buoyancy is smaller and 
the shear instability is favored. The resulting coefficient becomes ( Talon fc Zahnl . fl997t) 



-^shcar * '^*crit ^2 AT % 1 \^^) 



K+D h 



"■il.li ^T,ad = (Vad-V) and Nl = JLyV». (40) 



The horizontal turbulence favors the diffusion by the shear, reducing the inhibiting effect of the /i-gradient. 

There are other possible instabilities, which may play a role in different i ally r otating stars (baroclinic instabilities). 
The Goldreich-Schubert-Fricke (GSF) instability (|Goldreich fc Schubertl Il967ft which says that instability arises if 
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the T gradient, with account for thermal and viscous diffusivities, is insufficient to compensate for the growth of the 
centrifugal force during an arbitrary small displacement. If there is a significant /i-gradient, the instability is killed, 
thus it is generally not important. Ano t her p o ssible instability is th e so-called Axisymmetric-BaroClinic-Diffusivc 
(ABCD) instability dKnobloch fc Spruiti Il983t iTalon fe Zahil Il997h . If a fluid element moves over an arbitrary 
displacement from point A to C, it is brought back to C due to the conservation of angular momentum. However, in 
C the fluid eddy may be hotter than the surroundings (since the lines of constant f2 and T do not necessarily coincide 
in a baroclinic star). Due to thermal diffusivity, the fluid element loses some energy on its way. Thus, when back to 
A, the eddy is cooler and it overpasses its equilibrium position and the oscillation amplitude grows. The viscosity due 
to the horizontal turbulence tends to suppress this instability. On the whole, one must say that the real importance 
of the GSF and ABCD instabilities is still uncertain. 



V. ROTATION AND MAGNETIC FIELDS 

All the way from star formation to the endpoints, rotation and magnetic fields may be interacting. From the initial 
molecular clouds to the present Sun, the specific angul ar momentum dec reases by a factor ~ 10 6 . During the pre-MS 
phase, disk locking is a major effect for this reduction ( Hartmaiml Il998h . Fields of ~ 1 kG are sufficient for coupling 



a solar type star with an extended accretion disk. The contracting star is bound to its disk and it keeps the same 
angular velocity during contraction, thus losing lots of angular momentum. In solar type stars, the magnetic field 
creates a strong coupling between the star and the sol ar wind, which leads to further losse s of an gular momentum 
during the end of the pre-MS phase and the MS phase ( Kawaleii Il988t [Krishnamurthi et al. I. Il997t) . 



Spectropolarimetric surveys have obtained evid ences for the presence of magnetic field at the surface of OB stars 
(see e.g. the recent review by IWalder et all [20TTI and references therein). The origin of these magnetic fields is still 
unknown. It might be fossil fields (e.g. the spe ctral characteristics of Of?p stars are indicative of organized magnetic 
fields, most likely of a fossil origin according to I Wade et al l. l2010bl) . or fields produced through a dynamo mechanism. 
Recent simulations bv lCantiello et al I (120101) of subsurface convective zone in massive stars show dynamo-generated 



magnetic fields of the order of one kG. According to these authors, these magnetic fields might reach the surface of 
OB stars. 

A big question is whether there is a dynamo working in internal radiative zones. This could have far reaching 
consequences concerning the mixing of the elements and the loss of angular momentum. This question is also essential 
regarding the rotation periods of pulsars and the origin of GRBs (Sect. IVIII.DI and lXI.Bl) . It may also be important in 
order to explain the flat profile of the internal rotation of the Sun (jEggenberger et all 120051). Note that an a l ternat ive 



process has been invoked to explain this feature, the action of internal gravity waves ( Charbonnel fe Talonl l2005h 



A. Some properties of a dynamo in a radiative zone 

We examin e some properties of a dynamo in a radiative regi on. A particular case is t he Tayler-Spruit dynamo 
dSpruitl. 120021) and we are using many equa tions of this dynamo ( Maeder fc Mevnell . l2005al ). Some numerical simula- 



tions |Braithwaitel 120061 : IZahn et"aL f. I2007D confirm the existence of a magnetic instability, however the existence of 



the dynamo is still controversial. 

The two quantities we are mainly interested in for stellar evolution are the magnetic viscosity v, which expresses 
the mechanical coupling due to the magnetic field B, and the magnetic diffusivity r\ which expresses the transport by 
a magnetic instability and thus also the damping of the instability. Some general expressions for these two quantities 
may be derived. Let us first examine the energy conservation. The rate of magnetic energy production Wb per unit 
of time and volume is of the same order of magnitude as the rate W v of dissipation of the excess of rotational energy 
by the magnetic viscosity v. The differential motions in shellular rotation give velocity differences dv — rdQ. The 
amount of energy dissipated during a time dt for an element of matter dm in a volume dV is 

W v = \dm{dvf^j t = \ Q v{^ =^Oy with 9 = r|Vn|/n, (41) 

because the viscous time is dt = {dr) 2 /v. The magnetic energy density is u-q = B 2 /(8tt), it is produced within the 
characteristic growth time of the magnetic field cr^ 1 , thus the rate Wb of magnetic energy creation by unit of volume 
and time is 

B 2 1 

Wb = = -uj 2 A r 2 a B g , (42) 
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where we have used the expression of the Alfven frequency ua = r ^ ~p72 ■ Let us assume Wi/ = Wb, i.e. that 
the excess of energy in the differential rotation is converted to magnetic energy by unit of time. This gives for the 
viscosity coefficient of magnetic coupling 

This coefficient intervenes in the expression for the transport of angular momentum, in the Lagrangian form as given 
by Eq. (I27|) . Compared to the energy available for the solar dynamo, the amount of energy available from differential 
rotation is very limited. 

If due to an instability, some vertical motions with an amplitude l r /2 occur around an average position, the restoring 
buoyancy force produces vertical oscillations at the Brunt-Vaisala frequency N. The restoring oscillations will have 
an average density of kinetic energy % = / p v , wnere / is of the- order of unity. The magnetic held must overcome 
the buoyancy effect, which implies ub > otherwise the restoring force of gravity would counteract the magnetic 

2 

instability at the dynamical timescale. From this condition, one obtains l 2 < r 2 . If, / = A, we have the 



2/ N?- "> J 2 

condition (|Spruitl . 120021 ) 



lr<r^. (44) 

An unstable vertical displacement of size l r of the azimuthal held of lengthscale r and intensity B v also feeds a 
radial held component B r . The relative sizes of these two field components are defined by the induction equation, 
which gives the following scaling over the time St characteristic of the unstable displacement, B r rs SB « A B v St . 
For the maximum displacement l r given by Eq. (|44p . this gives 

W e a 7- < 45 > 

which provides an estimate of the ratio of the radial to azimuthal fields. 

Let us now consider the magnetic diffusivity 77, which damps the instability. The thermal diffusivity K produces 
heat losses from the unstable fluid elements and thus reduces the buoyancy forces opposed to the magnetic instability. 
Both effects have to be accounted for. If the radial scale of the vertical instability is small, the perturbation is quickly 
damped by the magnetic diffusivity rj (in cm 2 s _1 ). The radial amplitude must satisfy the relation, 

ll > — , (46) 

ctb 

where, as seen above, ctb is the characteristic frequency for the growth of the instability. The combination of the two 
limits (|4"4")) and (|4l)f gives for the case of marginal stability, 

" = t ■ < 47 » 

For given r\ and ctb, this provides the minimum value of wa, and thus of the magnetic field B, for the instability to 
occur. The instability is conhned within a domain, limited on the large side by the stable stratiheation (|44|) and on 
the small scales by magnetic diffusion (|4"r?)) . For the case of marginal stability, which is likely reached in evolution, 
this equation relates the magnetic diffusivity r\ and the Alfven frequency wa- 

The Brunt-Vaisala frequency N of a fluid element disp laced in a medium with account of both the magnetic and 
thermal diffusivities r\ and K is (M aeder fc MevneH |2004[ ). 



N 2 = ^~ A|, ad + A 2 , (48) 

K + L 

with iVf ad and N 2 defined in Eq. JUJ). The ratio -q/K of the magnetic to thermal diffusivities determines the heat 
losses. The factor of 2 is determined by the geometry of the instability , a factor of 2 applies to a thin slab, for a 
spherical element a factor of 6 is appropriate ( Maeder fc Mevnetl . l2005bl ). In the interior of a 15 M Q star, rj/K lies 



between 10 6 and 10 2 . Thus the T-stratification only has a little stabilizing effect, compared to the /^.-gradient. 
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Let us now examine the magnetic coupling, The momentum of force S by volume unity due to the magnetic field is 
obtained by writing the momentum of the Lorentz force F^. The current density j is given by the Maxwell equation 
i-£ j = v x B. Thus, one has 

S = f x F L = - f x (j x B) = —f x ( (V X B) X B) , (49) 

C 47T V / 

in modulus 5 « ± B r B v = J- B» = p r 2 (H) . (50) 

The units of S are g s~ 2 cm -1 , the same as for B 2 in the Gauss system. The kinematic viscosity v (in cm 2 s _1 ) for 
the vertical transport of angular momentum is 

V 1 „ dr 1 „ dr 1 „ c?lnr , . 

y _ J_ _ _ _ _ ~ — F f 511 

g g dv g rdfl g Cl d In fi ' 

where F is a force by surface unity, which also corresponds to a momentum of force by volume unity in g s~ 2 cm . F 
is applied horizontally to a slab of velocity v in a direction perpendicular to r. Considering only positive quantities, 
with q — |cilnf2/dlnr|, one has 

(52, 



pqVL NqQ 

Now, we can compare this expression for v to Eq. (|43[) and get 



UA^q 

= N , (53) 



which relates the growth rate of the magnetic field to its amplitude (through wa) 



B. The basic equations 

Introducing the expression ([531 of <jb in Eq. (|4"7|) . we get 



„ = r>ng(^) S . (54) 



Also, with the expression (|48|) , we can write for Cg 

, ,2 o 2 „2 
2 UJ A ll q 



-i+2 JV T, ad + J V 



(55) 



These equations are quite general. If the growth rate ctb of the instability is known, the two equations (|54p and ([55 
form a system of 2 equations with 2 unknowns 77 and wa ■ 



C. The case of Tayler-Spruit dynamo 

A pu r ely to roidal field B v (r, 1?), even very weak, in a stable stratified star is unstable on an Alfven ti mescale 1/oja 
(lTavled . [l97l . This is the first magnetic instability to appear. It is non-axisymmetric of type m = 1 ([Spruitj . Il999l 
I2002D . it occurs under a wide range of conditions and is characterized by a low threshold and a short growth time. For 
an azimuthal field consisting of concentric loops around the rotation axis, the instability appears as low-azimuthal 
order displacements of rings (Fig. [S]). Due to the magnetic pressure, the magnetic loops move apart from the rotation 

axis, like a disordered heap of tires, as de scribed by Spru it. The pressure is released by sideways motions. In a 

rotating star, the instability is also present ( Pitts fc Tavlerl Il985[ ). however the growth rate ctb of the instability is, if 



(56) 
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inste ad of the Alfven frequency u>a, because the growth rate of the instability is reduced by the Coriolis force ( Spruitl . 



2002). One usually has the following ordering of the different frequencies, N ^> >• cja- In the Sun, one has 
TV « 10~ 3 s _1 , f2 = 3 X 10~ 6 s _1 and a field of 1 kG would give an Alfven frequency as low as wa — 4 x 10~ 9 s _1 . 

The fre quency ctb being known, t he system of Eqs. (|54l) and (|55[) is defined and the 2 unknown rj and u>a can be 
obtained (jMaeder fe Mevnell l2005bh . When iV M dominates, which occurs around the stellar core, one has 

WA q IL an d ri = r 2 Qq i f—) 6 . (57) 



This shows that the mixing of the elements decreases strongly for large /j gradients and grows fast for large differential 
rotation. The ratio loa/Q has to be equal or larger than the mi nimum value defined by (J4TJ) . This leads to a condition 
on the minimum differential rotation for the dynamo to work ( Spruitl |2002T> , 



q > 



N\ 7/i / 77 \i/4 



When N 2 is large, as for example when there is a significant /x gradient, the differential rotation necessary for the 
dynamo to operate must also be large. 

Fig. [7] shows the differences of the internal f2-profiles during the evolution of a 20 Mq star with and without 
magnetic field created by the Tayler-Spruit dynamo. Without magnetic field, the star has a significant differential 
rotation, while is almost constant when a magnetic field created by the dynamo is present. It is not perfectly 
constant, otherwise there would be no dynamo. In fact, the rotation rapidly adjusts itself to the minimum differential 
rotation necessary to sustain the dynamo. 



VI. RECENT GRIDS OF MODELS OF MASSIVE ROTATING STARS 

In the recent years, many grids of massive rotating stellar models have been computed. A non-exhaustive list is 
given in Table HI These models have been computed accounting for the effects of rotation in the frame of the theory 
ex posed in the p revious sections. We have indicated those models which were computed with the the dynamo theory 
of ISpruitl (1999, 2002). Although the physics of rotation (and magnetic field when accounted for) is the same in all 
those models, the numerical implementations differ. Let us just mention here a few relevant differences. Some codes 
treat the transport of the angular momentum through a diffusion equation, not accounting for the advective nature 
of this equation. As explained above this may produce error not only on the amplitude of the effect but also on 
its sign! However, this statement should be tempered by the fact that for most of the MS phase, the region where 
the transport of angular momentum is transported efficiently, i.e. in the Gratton-Opik cell, angular momentum is 
transported outwards as would do a diffusion mechanism. 2) Rotating models contain uncertain parameters, at least 
two. One is in the expression of D s hear and account for the uncertainties pertaining our knowledge of the critical 
Richardson number (see Section IV. D). The other is in the expression of £>h, whose amplitude remains difficult to 
obtain from first principles. This coefficient has to be much greater than the one of D s hear in order to allow "shcllular 
rotation" to set in. Some models listed in Table U contain additional parameters, as a / M , a factor multiplying the 
mean molecular weight gradient. A small value of this parameter decreases the effects of the /i— barrier and facilitates 
mixing. To constrain the values of these uncertain parameters, some guidelines deduced from laboratory experiments 
or from multi-D dimensional simulations are used. Some authors also calibrate these uncertain parameters in order 
to achieve mean observed surface enrichments during the MS phase of B-type stars with models having an average 
rotational velocity. 

A few authors have built 2D models for rotating stars ( Deupreel .[T990; Es pinosa Lara fc Rieutordl 20071 : iRoxburghl 



2006). For instance, ZAMS models with arbitrary rotation laws have been computed bv lDeupreej "( 200lh . The variation 



with the latitude o f the effective temperature for rotationally deformed stars from a 2D stellar s tructure is given by 
lLovekin et al. (2006). Non radial oscillations are computed for 2D rotating massive mod els by lLovekin fc Deupred 



( 200S). Effect s of ro tation on the pulsation frequencies has been investigated recently bv lDeupree fc Beslinl ( 2*0101 ); 



lLovekin et all (|2009h 



VII. THE ELEMENT ABUNDANCES IN MASSIVE STARS 

Both mass loss and rotational mixing affect the evolution of the chemical abundances at the stellar surface of 
massive stars. At solar abundances, mass loss dominates in stars initially more massive than about 30 M Q , while 
rotation dominates below this mass limit. Let us first examine the effect of mass loss in the most massive stars. 
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A. Model effects of mass loss on abundances 

Fig. [8] shows the internal evolution of a star with an initial 60 M Q . The removal of the outer layers reveals the 
internal composition. There are 5 typical domains of composition, which result from a progression in the exposition 
of nuclear products. 

- 1. The initial abundances. Mass loss does in general not change the surface composition in OB stars during the MS 
phase, except for stars above 60 M Q . 

- 2. Intermediate abundances. This stage is characterized by partial CNO processing with possible dilution effects. 
There are N enrichments, enhancements of the 13 C/ 12 C ratio, C depletions and modest O depletions. This is typical 
of blue and red supergiants. 

- 3. CNO equilibrium with H present. CNO equilibrium is reached before H exhaustion. The C/N and O/N ratios 
are reduced by two orders of a magnitude with respect to cosmic abundances. This is the case of WR stars (bare 
cores) of types WNL. 

- 4. CNO equilibrium with H absent. The He mass fraction is 98%, the CNO ratios are the same as before. This 
stage corresponds to WR stars of type WNE. The abundances are model independent, being determined mainly by 
nuclear cross-sections. 

- 5. Partial He burning. The products of He burning, i.e. 12 C, 16 0, 22 Ne are visible. The changes are abrupt 
(rotation make them smoother). The abundances depend strongly on the models (mass loss, mixing, etc.). This stage 
corresponds to WC stars and to WO stars for O/C ratios > 1. 

Of course, not all stars go through this whole sequence: the smaller the initial mass and/or metallicity, the shorter 
the path. Rotation may accelerate and make smoother the evolution of surface composition, due to the contribution 
of internal mixing. 



B. The effects of rotation on the surface abundances 

In the mass range below about 30 M©, mass loss has little effect during the MS phase. The changes of abundances, 
if any, are due to mixing. Fig. [5] shows the run of the diffusion by shear turbulence, -D s hcar, and by the net effect of 
both meridional currents and horizontal shear turbulence, -D e ff, at three different stages during the MS phase of a 
rotating 20 M Q model. The main mixing effect is the diffusion by shear turbulence, which results from the internal Q 
gradients built during evolution. To a smaller extent, meridional circulation makes some transport, however mainly 
of angular momentum. 

Mixing brings to the surface the products of CNO burning: mainly 14 N and 13 C enrichments, 12 C is depleted with 
limited 4 He enrichment and 16 O depletion. Fig. [TO] shows the predicted variations of log(N/H) during MS evolution as 



a function of the initial masses ([Mevnet fc Maederl . 120001) . (N/H) is here the abundance ratio of N and H in numbers. 



Without rotational mixing, there would be no enrichment until the red supergiant stage, but rotation produces an 
increase (depending on velocity v) of N/H during the MS phase. The N excesses also depend on the ages t. The 
increase is modest during the first third of the MS phase, because the elements need some time to reach the surface, 
then it is more rapid. The N enrichments are larger for larger masses M . 

Fig. HU shows, for given M at Z = 0.02, the evolution of the N excesses with age for different initial velocities. It 
shows that there is no single relation between vsini and the N/H excesses for a mixture of ages. The same is true 
for a mixture of different masses. Binarity may also affect the N and He enrichments due to tidal mixing and mass 
transfer. A binary star with low rotation may have a high N/H due to tidal mixing or due to the transfer of the 
enriched envelope of a red giant. At the opposite, a binary star may also have a high vsini and no N/H excess, in 
the case of the accretion of an unevolved envelope bringing a lot of angular momentum. On the whole, the N excess 
is a multivariate function 

Alog(N/H) = f(M, t, vsini, multiplicity, Z) . (59) 



Models with lower initial metallicities Z have higher N enrichments for given M and v ( Maeder fc Mevnetl . l2001at : 



iMevnet et a"ill2006f ). The reason is the higher fi-gradients, resulting from the absence of the Gratton-Opik circulation 
cell (Sect. IIV.B|) . The excesses become very strong at metallicities Z as low as 10~ 8 . 



C. The observed N/H excesses 

The amplitudes of the N enrichments at the end of the MS phase in massive stars form a reference point telling 
us the importance of mixing. The data at different Z (mean and largest values) are summarized in Table [H] from a 
number of sources (see references below the Table). In the lowest mass range considered (6.6-8.2 M Q ), small excesses 
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of He/H are still present (jLvubimkov et al. 
We see the following facts: 



. 120041 ). they are larger in the Small Magellanic Cloud (SMC, Z « 0.004). 



• On the average, the N enrichments are larger for larger masses. 

• The N enrichments are larger at lower Z. 

• Away from the ZAMS, but still in the Main Sequence, the He and N enrichments are larger and they are even 
larger in the supergiant stages. These various features are quite consistent with the predicted properties of 
rotational mixing as can be seen in Fig. [T2] 

• Correl ations between N or He excesses and the observed v sin i have been found in the upper part of the MS 
band (|Lvubimkov et all 12004ft , in agreement with model predictions. 



To find a correlation for a multivariate function like N/H with some parameter like vsini, it is necessary to limit 
as much as possible the range of the other involved parameters. Otherwise, the conclusions may be erroneous. From 
Fig. 1131 we note that data samples limited in masses and ages support a N enrichment depending on rotational 
velocities. Stars be yond the end of th e MS phase do not obey to such a relation, because their velocities converge 
toward low values ( Mevnet fc MaedeH [2000) . A fraction, which we estimate to be less than ~ 20 % of the stars, 
may esc ape to the relation as a result of binary evolution or due to some other mechanisms as for instance magnetic 
braking ( Mevnet et al.l. 2011 ). This l atter mechanism might explain some magnetic stars showing large N-excesses 
and low v sin i dHenrichs etall l2003al lbh . 

To obtain more tight constraint s on the models, very careful and accurate abundance analysis are required. Very 
interestingly iPrzvbilla et al] ( 2010h showed that observed ratios of N/O and N/C follows a well defined relation in the 
N/O versus N/C plane, which does not depend on the models but just on the CNO nuclear reactions properties. In 
the future, this characteristic can be used to test the surface abundances obtained from spectral analysis before these 
abundances are used to constrain t he stellar models. 

According to lHunter et al.l ( 20071 ) main-sequence binary objects have close to baseline nitrogen surface abundances. 
These systems thus do not present apparent signs of extra-mixing. In contrast several evolved binary objects have 
high nitrogen enhancements. These abundances are similar to those observed in apparently single stars. Thus it 
appears difficult to discriminate among the possible causes of the enrichments in binaries, i.e. between extra- mixing 
operating in s i ngle s tars and mass-transfer events in close binary systems. The same result has been obtained by 
iTrundle et ail f)2007h . 



D. The observed He excesses 

Surface enrichments in helium have been observed with the following main trends: 



He in O-type stars in the SMC: In the SMC, for 31 O-type stars, iMokiem et~aH (|2006h find values of 



y 



nn e /{riYi + riHe) between 0.09 and 0.24, where hhc is the density number of helium and «h of hydrogen . Note 
that riHe/ (^h+^Hg) = Y/ (Y+AX), where Y and X are respectively the mass fraction of helium and of hydrogen. 
Setting Y ~ 1 — X, one obtains values for Y = 4y/(l + 3y) between 0.28 and 0.56 (here 0.28 would correspond 
to y=0.09 i.e. to the initial helium mass fraction). These authors conclude that while rotation can qualitatively 
account for such enrichments, the observed enrichments are in many cases much stronger than those predicted 
by the models. 

He in O-type stars in the Large Ma gellanic Cloud ( LMC ): In the LMC, for 28 O-type stars, y values between 
about 0.09 and 0.28 are obtained by Mo kiem et al.l (|2007t) . i.e. helium mass fractions between 0.28 and 0.61. 



He in OB-type stars in the Milky Way (MW): lHuang fc Gied (|2006bl ) determine He abundances for OB galactic 



stars. In their high mass range (8.5 M Q < M < 16 M Q ), the He enrichment progresses through the main sequence 
and is greater among the faster rotators 1 . On average He abundance increases of 23% ±13 % between ZAMS 
and T AMS. These authors also obtain that He enrichments are higher for higher v sini values. iLvubimkov et all 
(2004) :ind an increase during the MS phase in He abundance of 26% for stars in the mass range 4-11 M Q and 



67% for more massive stars in the range 12-19 M 



1 These authors also found many helium peculiar stars (He- weak and He-strong). These stars were not used to study the process of 
He-enrichment. 
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From a theoretical point of view, one expects that the N-enrichment occurs very rapidly at the surface well before 
any surface helium enrichment. This is of course due to the fact that very rapidly the nitrogen abundance increases 
at the centre creating thus a strong chemical gradient between the core and the envelope. Since diffusive velocity 
is greater when the gradient of abundance is greater, the presence of such a strong gradient favors a rapid mixing. 
The gradient of helium in stellar interiors is built up on much longer timescales than the gradient in nitrogen. As a 
consequence, the diffusion of this element is expected to be much less rapid than that of nitrogen. Thus any strong 
He-enhancement should be accompanied by a very strong nitrogen enhancement. At the moment, the observations of 
both He- and N-enrichments in the same stars are too scarce to check this trend. 



VIII. EVOLUTION OF ROTATION AT THE SURFACE AND IN THE INTERIOR 

A. Initial conditions for rotating models 

Ideally one would like to know the shape of the distribution of rotational velocities on the ZAMS for the whole 
mass and metallicity ranges spanned by the stars during the whole cosmic history. This is of course not an easy task 
and many difficulties have to be overcomed. Let us just mention the main ones here: 

• First, most observations concern non-ZAMS stars, whose surface velocity is not representative of the initial 
stellar velocity. In order to obtain ZAMS velocities some theoretical guidelines have to be used. For instance, 
using stellar models, it is possible to deduce which are the velocities needed on the ZAMS in order to achieve a 
given mean velocity during the MS phase. 

• What is measured through Doppler broadening is vsmi and thus a sufficiently great number of stars need to 
be measured in order to correct in a statistical way for the effect of the inclination angle. 

• There are theoretical reasons for believing that beyond a velocity limit the line broadening due to rotation 
saturates, making this method inappropriate for measuring velocities near the critical limit. This is because 
when the velocity approaches the critical one, the equator, which has the higher linear velocities, darkens and 
thus contributes less to the line formation than the regions of the star near to the pole which have the smaller 
linear velocities. 

It is also worthwhile to underline the fact that even if we knew very well the distributions of the surface velocities, 
this would not help in determining the initial distribution of f2 inside the star. Hopefully asteroseismology wil l allow 
us to obtain such information, but presently the indications are still scarce (see however results presented in lAertsl . 
2008) 2 . Also interferometry, which allows us to measure the shape of fast rotating stars, can give some indications on 
how the interior rotate. For instance a ratio of the equatorial radius to the polar radius equal to 1.5 when the star 
rotates at the critical limit implies that the interior is not much deformed, in agreement with the Roche approximation. 

The knowledge of the interior distribution of f2 would be important for estimating the total angular momentum of 
the star on the ZAMS. At present time, stellar models make the simplifying hypothesis that the evolution starts with 
a flat profile of f2. The exact shape of the initial profile of f2 is not so important for the further evolution of f2 inside 
the star. It has been found that the initial distribution of the angular momentum ins ide the star is rapidly for gotten 
and erased by a rapid convergence of the profile of f2 towards an equilibrium profile ( Denissenkov et all Il999f) . This 



results from two counteractive effects: on one hand meridional currents build fl gradient (which tends to increase Q 
in the inner layers) and on the other hand shear turbulence tends to decrease the fi gradient. 

The mean initial angular momentum of a model star on the ZAMS, Ani) is estimated as the value needed for the 
star to rotate at the observed average rotational velocity during the MS phase. The inertia momentum of a 10 M Q 
star on the ZAMS is of the order of 10 56 g cm 2 . Considering a radius R ~ 4i? Q and a surface equatorial velocity, v eq , 
of 200 km s _1 (see Fig. rjj, one obtains for f2 = v eq /R = 0.00007 s _1 , and a total angular momentum content of the 
order of 10 52 g cm 2 s _1 assuming the star rotates as a solid body on the ZAMS. These are the typical initial angular 
momentum which are considered in stellar models. 



2 See Table 1 in this reference. Two stars present a high contrast between the angular velocity of the core and that of the envelope, 
between a factor 4 and 5. One star does not show any contrast. All the three stars are slow rotators and are B2-B3 stars with class 
luminosity between V and III. 
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B. Evolution of the surface velocities 



Once the initial conditions have been chosen, models predict the evolution of the surface velocity as a function of 
time. The evolution of v in the HR diagram is shown in Fig. Q3] (v ini — 300 km s _1 on the zero-age main sequence). 
Lines of constant v are indicated. On the MS, we notice the decrease of the surface velocity as a function of time. 
This trend is confirmed by the observations of iHuang fe Giesl ( 2006bl ) who show that all OB stars of their sample 



experience a spin-down during the MS phase. A few relatively fast rotators are found near the end of the MS phase. 
According to these authors, these stars may be spun up by a short contraction phase or by mass transfer in a close 
binary. In Fig. [14j we see also that the decrease of the surface velocity is much faster for the most massive stars than 
for stars with M < 15 M©, as a consequence of stronger mass losses in th e higher mass range . This difference remains 
also present in the domain of B-supergiants. Interestingly, according to iHunter et al.l (120081) the r e is so me evidences 



that the most massive objects rotate slower than their less massive counterparts. Dufton et al.l ( 20061 ) confirm that 



the mean rotational velocity of stars which have strong winds is lower than that of the lower mass stars. 

A few other points can be noted: the average v is lower for O-type stars than for the early B-type stars ( Slettebakl . 



Il970h . A gain this may be the consequence of the higher losses of mass and angular momentum in the most massive 
stars. Also, we r emark that the increase of v from O-stars to B-stars is larger for the stars of luminosity class IV 
than for class V (iFukudal . Il982h . This is consistent with the models, which show (cf. Fig. [T4|) that the differences 



of v beween O- and B-type stars are much larger at the end of the MS phase. Another fact in the observed data 
is the strong decrease of v for the massive supergiants of OB-types. Also during the crossing of the HR diagram, 
the rotational velocities decrease fast, to become very small, i.e. of t he order of a few km s _1 , in the red supergiant 
phase. This is predicted by all stellar model s (cf. also Langer . 19981) due to the growth of the stellar radii. This is 



also confirmed by the recent observations of iMokiem et al. I (|2006l) and lHunter et al. (2008). These last authors find 



from a sample of 400 O- and early B-type stars in the Magellanic Clouds that supergiants are the slowest rotators in 
the sample, typically having rotational velocities less than 80 km s -1 . 

Note that the models shown in Fig. [T?] are models without magnetic fields. When magnetic fields are accounted for 
as in the Tayler-Spruit dynamo, the evolution of the surface velocity is different. For instance, a 15 M© model with 
Wini=300 km s _1 on the ZAMS and computed with magnetic fields, keeps a nearly constant surface velocity around 
300 km s _1 all along the MS phase . Without ma gnetic fields, the velocity decreases during the MS phase. The 
average velocity is around 240 km s _1 (see Fig. 3 in lMaeder fe Mevnett |200~4| ). 

Fig. [15] clearly illustrates the very different evolution of the rotational velocities of a 60 M© at various metallicities. 
At low Z like in the models at Z = 10~ 5 , the growth of ^- is possible because of the very small mass loss by 
stellar winds. In view of these results, it is likely that at very low Z a large fraction of the massive stars reach their 
break-up velocities (see right panel of Fig. |2"4"|) . This question is of high importance, because if the massive stars 
reach their break-up velocity, most of their evolutionary and structural properties will be affected. For example, they 
could also lose a lot of mass and produce some Wolf-Rayet (WR) stars. They would have a relatively small remaining 
mass at the time of the superno va explosion, l ike th eir counterparts at solar composition. Interestingly, the trend 
illustrated above is supported bv lHunter et al.l (l2008f) who obtain that SMC metallicity stars rotate on average faster 
than galactic ones (mainly field objects) and by iMartavan et all ( 20071 ). who find that, for B and Be stars, the lower 
the metallicity, the higher the rotational velocities. 



C. The Be stars and the evolution towards the critical velocity 

Be stars are emission line stars. Emission originates in a circumstellar outflowing disk. How do these disks form ? 
How long are their lifetimes ? Are they intermittent ? Are they Keplerian ? Many of these questions are still subject 
of a lively debate. A point however which seems well accepted is the fact that the origin of a disk might be connected 
to the fast rotation of the star. This view is support ed by the fact that objects with Be phenomena are the fastest 
rotators in the sample studie d bv lHunter et al.l (2008) (400 OB stars in the MCs). This trend has also been found by 



IMartavan et al.1 ( 2007 . 120061) who obtain that Be stars rota te faster than B stars whatever the metallicity. 



Maeder et al.l ( 1999 ) and lWisniewski fc Biorkmanl (|2006l) find that the fraction of Be stars with respect to the total 
number of B and Be stars in clusters with ages (in years) between 7.0 and 7.4 (in logarithm) increases when the 
metallicity decreases. This fraction passes from about 10% at solar metallicity to about 35% at the SMC metallicity. 
These significant proportions imply that Be stars may drastically affect the mean v sin i obtained for a given population 
of B-type stars. Very interestingly there appears to be a correlation between the frequency of Be stars and that of 
red supergiants. 

These observations indicate that metallicity plays a role in the Be phenomenon and provides hints on the way 
surface velocity may evolve differently for stars of different initial metallicities. The correlation of Be star populations 
with those of red supergiants can be seen as an indication that fast rotation not only favors the formation of Be stars 
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but also that of red supergiants (see Sect. IX.A[) . 

The evolution of surface velocities during the Main Sequence lifetime results from an interplay between meridional 
circulation (bringing angular momentum to the surface) and mass loss by stellar winds (removing it). The dependence 
on metallicity of these two mechanisms plays a key role in determining for each me tallicity, a limit i ng rang e of initial 
masses (spectral types) for stars able to reach or at least approach the critical limit. lEkstrom et al.l (|2008bl) computed 
models (with no account for the Tayler-Spruit dynamo) for initial masses between 3 and 60 M© and considered, for 
each initial mass, models with f2/f2 cr it equal to 0.1, 0.3, 0.5, 0.7, 0.8, 0.9 and 0.99. The following conclusions have 
been obtained based on these models: 

1. Meridional currents bring angular momentum inwards only at the very beginning of the core H-burning phase. 
During most of the MS phase, meridional currents transport angular momentum from the inner regions to the 
outer layers. As long as the mass loss rates are not too important, this brings the v/v C nt ratio close to 1 during 
the MS phase for stars having on the ZAMS f2/r2 cr ;t >^ 0.8. 

2. Since the velocity of the meridional currents in the outer layers scales with the inverse of the density, the process 
becomes more efficient for stars of higher initial mass and/or higher initial metallicity. 

3. At high metallicity however, mass loss becomes more and more important and can prevent the stars from 
reaching the break-up limit. 



4. In stellar clusters, supposing an initial distribution of the rotational velocities as given bv lHuang fc Giesl ( 2006ah . 



one expects that the fraction of stars having u/w C rit > 0.8 becomes maximum for ages between 20-32 Myr at 
Z = 0.020. This range of ages shifts to older ages at higher metallicities and to younger ages at lower metallicities. 
Such a shift results from the dependence on metallicity and mass of the intensity of the mass loss rates and of 
the velocities of the meridional currents. If we look at stars with v/v CI - lt > 0.7, the maximum number will be 
found in clusters of ages around 10 Myr and below, and this at all non-zero metallicities. 

5. Be stars might be the natural outcome of stars with an initial rotational velocity in the upper tail of the initial 
velocity distribution. 

6. Depending on when the critical velocity is reached, one expects more or less high surface enrichments (see Fig. 
I12[) . If the limit is reached very early during the MS phase, no enrichment is expected, while if this limit is 
reached at the end of the MS phase, high N/C and N/O ratios are expected. However, the enrichment also 
depends on the initial mass. For stars originating from small initial mass stars (typically from 3 M Q at standard 
metallicity) one expects no or small surface enrichments, at whatever time the critical velocity is reached. 

7. To reproduce the higher fractions of Be stars at low metallicity, we must assume that in metal poor regions, a 
higher number of stars are born with high values of f2/f2 cr it- 



D. Evolution of the interior angular momentum content 

Figu re IT6| shows the evo lution of fi inside the 25 M Q model from the ZAMS until the end of the core Si-burning 



phase ( Hirschi et all I2004T ) . The evolution of fl results from many different processes: convection enforces solid body 



rotation, contraction and expansion respectively increases and decreases f2 as a consequence of local conservation of 
the angular momentum, shear (dynamical and secular) erodes f2-gradients while meridional circulation may erode or 
build them up and finally mass loss may remove angular momentum from the surface. If during the core H-burning 
phase, all these processes may be important, from the end of the MS phase onwards, the evolution of O is mainly 
determined by convection, the local conservation of the angular momentum and, for the most massive stars by mass 
loss. 

During the MS phase, ft decreases in the whole star. When the star becomes a red supergiant (RSG), il at the 
surface decreases significantly due to the expansion of the outer layers. Note that the envelope is gradually lost by 
winds in the 25 M a model. In the centre, fl significantly increases when the core contracts and then the 51 profile 
flattens due to convection. f2 reaches values of the order of 1 at the end of Si-burning. It never reaches the local 
break-up angular velocity limit, 51 c , although, when local conservation holds, fl/fl c oc r -1 / 2 . 

Figure [TTl shows the evolution of the specific angular momentum, j — 2/3 51r 2 , in the central region of a 25 M 
stellar model. The specific angular momentum remains constant under the effect of pure contraction or expansion, 
but varies when transport mechanisms are active. One sees that the transport processes remove angular momentum 
from the central regions. Most of the removal occurs during the core H-burning phase. Still some decrease occurs 
during the core He-burning phase, then the evolution is mostly governed by convection, which transports the angular 
momentum from the inner part of a convective zone to the outer part. This produces the teeth seen in Fig. 1171 The 
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angular momentum of the central regions of the star at the end of Si-burning is essentially the same as at the end 
of He- burning (by end of He-burning, we mean the time when the central helium mass fraction becomes less than 
10 -5 ). As a numerical example, we calculated, for the 25 M Q model, the angular momentum of its 3 M Q remnant. We 
obtained £, em = 2.15 10 50 gcm 2 s _1 at the end of He-burning and £ rom = 1-63 1 50 gcm 2 s _1 at the end of Si-burning. 
This corresponds to a loss of only 24%. In comparison, the angular momentum is decreased by a factor ^5 between 
the ZAMS and the end of He-burning. This shows the importance of correctly treating the transport of angular 
momentum during the Main Sequence phase. It means also that, in the frame of the present models, we can obtain 
a first order of magnitude estimate of the pre-supernova angular momentum of the core at the end of He-burning. 

The p resent model s have difficultie s in accounting for the ob served rotation rates of young pulsars and white dwarfs 
(see e.g. iHeger et all 120051: iKawaieR 119881: ISuiis" et all 2008). More precisely, they predict too fast rotation of the 
stellar cores in the advanced phases. This may be due to different causes that could lead to additional angular 
momentum loss from the central regions at different evolutionary phases: 

• already during the nuclear lifetime; 

• at the time of the supernova explosion in the case of neutron stars, or during the Thermal Pulse Asymptotic 
Giant Branch (TP-AGB) phase at the time of the superwind episode in the case of white dwarfs; 

• during the early evolution of the new born neutron star or white dwarf. 

During the MS phase, angular momentum can be efficiently extracted from the core in case the star rotates as a solid 
body. Such a situation is realized in numer i cal sim u lations where the T ayler-Spruit dynamo is accounted for as can be 
seen in Fig. [7] (see e.g. iMaeder fc Mevnetl . l2005bl ). IHeger et al.l (2005) find that magnetic torques decrease the final 



rotation rate of the collapsing iron core by about a factor 30-50 when compared with the nonmagnetic counterparts, 
allowing to obtain rotation periods of young pulsars which better fit the observations. 



IX. EVOLUTIONARY TRACKS AND LIFETIMES 

A. The Main-Sequence phase 

Figure [T5] shows the evolutionary tracks of non-rotating and rotating stellar models for initial masses between 9 and 
120 M© at Z=0.020. For the rotating stellar models, the initial velocity is v- m i = 300 km s _1 . There is little difference 
between tracks with w ini = 200 or 300 km s _1 . If the effects behaved like vf ni , there would be larger differences. 
The present saturation effect occurs because outward transport of angular momentum by shears are larger when 
rotation is larger, also a larger rotation produces more mass loss, which makes a larger reduction of rotation during 
the evolution. Let us note however that for some surface abundance ratios as N/C or N/O, the increase from Vi n i = 
200 to 300 km s _1 produces significant changes. Thus, the similarity of the evolutionary tracks does not necessarily 
imply the similarity of the surface abundances for these elements. 

On and near the ZAMS, rotation produces a small shift of the tracks towards lower lumin osities and T e g. This 
effect is due to both atmospheric distorsions and to the lowering of the effective gravity (see e.g. lCollins fe Sonnebornl . 
ll977tlKippenhahn fc Thomasl . ll970l : lMaeder fc Pevtremannl . fi 9 70() . At this stage the star is still nearly homogeneous. 
When evolution proceeds, the tracks with rotation become more luminous than for no rotation. This results from 
essentially two effects. On one side, rotational mixing brings fresh H-fuel into the convective core, slowing down its 
decrease in mass during the MS. For a given value of the central H-content, the mass of the convectiv e core in the 
rotat i ng model is the r efore larger than in t he non-rotating one and thus the stellar luminosity is higher ( Heger et all 
120001: lMaede"rl . Il987bl: iTalon fc Zahnl Il997l) . 

Rotation reduces the MS width in the high mass range (M ~ 40 M©). Let us recall that when the mass increases, 
the ratio of the diffusion timescale for the chemical elements to the MS lifetime decreases. This can be shown in the 
following way: let us call Tdiff the diffusion timescale which is proportional to R 2 /D where R is the radius and D 
the diffusion coefficient. For shear mixing, D is proportional to K — AacT 3 / (3Kp 2 cp) , the thermal diffusivity. Using 
homology relations, it can be shown that Tdiff varies as 1/M 1 - 8 . Now the MS lifetime, tms, for stars with initial 
masses greater than 15 M© varies as 1/M 8 , therefore the ratio Tdis/rMS varies as M~ 1A . As a consequence, starting 
with the same «i n i on the ZAMS, massive stars will be more mixed than low mass stars at an identical stage of their 
evolution. This reduces the MS width since greater chemical homogeneity makes the star bluer. Moreover, due to 
both rotational mixing and mass loss, their surface will be rapidly enriched in H-burning products. These stars will 
therefore enter the Wolf-Rayet phase while they are still burning their hydrogen in their core. This again reduces the 
MS width. For i nitial masses betw een 9 and 25 M Q , the MS shape can be widened by rotation mimicking the effect 
of overshooting ()Talon et all Il997t ). 
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B. The post- Main-Sequence evolution 

The post-MS evolution of the most massive stars (M > 40 M Q ) which become WR stars will be discussed in Sect. 
IX. Bl below. We mention here some points of general interest. For low or moderate rotation, the convective core 
shrinks as usual during MS evolution, while for high masses (M ~ 40 M Q ) and large initial rotations (fp— > 0.5), the 
convective core grows in mass during evolution. This latter situation occurs in the fast rotating 60 M Q model shown 
on Fig. 1181 These behaviours, i.e. reduction or growth of the core, determine whether the st ar will fo l low res pectively 
the usual redwards MS tracks in the HR diag ram, or whether it w ill bifurcate to the blue (jMaederL Il987bh towards 
the classical tracks of homogeneous evolution ( Schwarzschildl . [19581 ) . 



The stars with initial masses between 15 and 25 M become red supergiants (RSG). Rotation does not change 
qualitatively this behaviour but accelerates the redwards evolution, especially for the 15 and 20 M Q models. As a 
numerical example, for an initial v- ln i = 300 km s _1 , the model stars burn all their helium as red supergiants at T e s 
below 4000 K, while the non-rotating models spend a significant part of the He-burning phase in the blue part of the 
HR diagram: for the non-rotating 15 and 20 M e models, respectively 25 and 20% of the total He-burning lifetime is 
spent at logT c ff > 4.0. Note that rotation is not the only factor able to affect the blue-red supergiant lifetimes. Mass 
loss, convection, overshooting, semiconvection, are, among others, all important factors in that respect. The above 
behaviour of the rotating models results mainly from the enhancement of the mass loss rates. This effect prevents 
the formation of a big intermed i ate co nvective zone and therefore favours a rapid evolution toward the RSG phase 
(jMaedeii 1 1 9 8 it Istothers fc Chinl . Il979h . Let us note that the dispersion of the initial rotational velocities produces a 
mixing of the above behaviours. As we shall see below (see Sect. IX.Aft . a redwards evolution is also favored at low 
metallicity. But in that case, this is main ly due to rotational mixing and not to a mass loss effect. 

The problem of the mass discrepancy ( Herrero et all li~992l ). i.e. of the difference obtained between spectroscopic 



and evolutionary masses, has been significantly reduced thanks to improvements brought to stellar atmosphere models. 
However still some discrepancies are reported in general linked with strong helium surface enrichments: 

• In SMC: lMokiem et al. 1 (120061) find a mild mass discrepancy for stars with spectroscopic masses inferior to about 



20 M Q , which correlates with the surface helium abundance. These authors find that the discre pancies are 
consis tent with the predictions of chemically homogeneous evolution. Most of the stars observed bv lHeap et al.l 
(2006) exhibit the mass discrepancy problem although no surface He enrichment. 



In LMC: iMokiem et al. I (120071) from the analysis of O-type stars in the LMC find that bright giants and su- 



pergiants do not show any mass discrepancy, regardless of the surface helium abundance. In contrast they find 
that the spectroscopically determined masses of the dwarfs and giants are systematically smaller than those 
derived from non-rotating evolutionary tracks. All dwarfs and giants having y > 0.11 (Y > 0.33) show this mass 
discrepancy. 

These mass discrepancies may arise from the use of stellar models with too small convective cores. Convective 
overshooting and rotation produce larger convective cores and thus can help in explaining this mass discrepancy. 

C. Lifetimes 

Generally we can say that the MS lifetime duration is affected by rotation at least through three effects: 

1. Rotation increases the quantity of hydrogen burnt in the core. This increases the MS lifetime. 

2. The hydrostatic effects of rotation make a star of a given initial mass to behave as a non-rotating star of a 
smaller initial mass. This tends to increase the MS lifetime. 

3. Rotation increases the helium abundance in the outer radiative envelope. This tends to make the star overlu- 
minous with respect to its non-rotating counterpart and thus to reduce the MS lifetime. 

Depending on the initial mass and metallicity, rotation can either increase or decrease the MS lifetimes with respect to 
the durations obtained from non-rotating stellar models. For Z = 0.020, the e ffect n umber 1 dominates and rotation 
increases the lifetimes. For instance, from the models by iMevnet &: Maederl ( 2000h . one can deduce a nearly linear 



relation between the relative enhancement of the MS lifetime, A% and v, where 

A%(w) = M«)-tH(0)]/tH(0). 

One obtains, 



0.0013 ■ v, 
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where v is in km s _1 . When the metallicity decreases, the effect number 3 tends to become the most important one. 
Typically at Z = 10 -5 , the MS lifetimes are decreased by about 4 -14% for the mass range between 3 and 60 M 
when w; n i increases from to 300 km s _1 ( Mevnet fe Maederl l2002al ). The rotating 2 M Q model of the same grid has 
a longer MS phase than its non-rotating counterparts, in contrast with what happens for higher initial mass stars. 
This is because, when the initial mass decreases, the hydrostatic effects (effect number 2) become more and more 
important. 

The He-burning lifetimes are less affected by rotation than the MS lifetimes. The changes are less than 10%. The 
ratios tHe/tn of the He to H-burning lifetimes are only slightly decreased by rotation and remain around 10-15%. 



X. MASSIVE STAR POPULATIONS AND THE CASE OF FIRST STARS 
A. The ratio of blue to red supergiants 



As noted bv lLanger fc Maederl (1995), the current models (without rotation) with Schwarzschild's criterion predict 



no red supergiants in the SMC. This is also seen in Fig. H~9l which illustrates for models of 20 M Q at Z = 0.004 the 
variations of the T e s as a f unction of the fractional lifetime in the He-burning phase for different rotation velocities 
(|Maeder fc Mevnetl . [2001bh . 7 or zero rotation, we see that the star only moves to the red supergiants at the very end 



of the He-burning phase, so that the blue to red supergiant ratio B/R is ~ 47. This is 50-100 times greater than the 
ob served ratio in the SMC clust er NGC 330 which l i es bet ween 0.5 and 0.8, according to the various sources discussed 
in lLanger fc Maede"rl ( 19951) and lEggenberger et al I 1)20021) . 



This disagreement implies that no reliable predictions can be made concerning the nature of the supernova pro- 
genitors in different environments, or the populations of supergiants in galaxies. The B/R ratio also constitutes an 
important and sensitive test for stellar evolution models, because it is very sensitive to mass loss, convection and mix- 
ing processes ( Langer fc Maederl . Il995l ). Thus, the problem of the blue to red supergiant ratio (B/R ratio) remains 



one of the most severe problems in stellar evolution. 

Looking at Fig. QjU we see that for average rotational velocities v during the MS, v — 152, 229 and 311 km s -1 , 
one has respectively B/R = 1.11, 0.43 and 0.28. Thus, the B/R ratios are much smaller for higher initial rotation 
velocities, as rotation favours the formation of red supergiants and reduce the lifetime in the blue. We notice in 
particular that for v — 200 km s _1 , we have a B/R ratio of about 0.6 well corresponding to the range of the observed 
values. 

The reason for t his va riation of the B /R ratio when the rotation velocity changes is discussed in details in 
IMaeder fc Mevnetl (l2001bD . Let us just mention here that when an extended intermediate convective zone is as- 



sociated to the H-burning shell, the star tends to stay in the blue part of the HR diagram Rotational mixing tends 
to reduce the extension of this intermediate convective zone and thus to accelerate the redwards evolution in the HR 
diagram. How rotation limits the convective zone associated to the H-burning shell can be understood looking at 
Fig. [20] The larger He-core in the rotating models means that a larger fraction of the total luminosity is made in 
the core (0.42 instead of 0.31 for the models of Fig. |2"0"1) . This means that the H-burning shell in the rotating model 
produces a smaller fraction of the total luminosity and this contributes to reduce the importance of the convective 
zone above the H-burning shell. Simultaneously, the higher He-content near the shell in the rotating case leads to 
a decrease of its H-content and opacity, and this also contributes to reduce the importance of the convective zone 
associated with the H-burning shell. 

Thus models with rotation well account for the long standing problem of the large numbers of red supergiants 
observed in low Z galaxies, while current models with mass loss were predicting no red supergiants. 



B. Current problems with Wolf-Rayet stars 

Wolf-Rayet stars are the bare cores of initially massive stars, whose H-rich envelope has been removed by strong 
stellar winds or through Roche lobe overflow in a close binary system (see e.g. the review bv lCrowtheitl2007ri . Wolf- 
Rayet stars play a very important role in Astrophysics, as signatures of star formation in galaxies and starbursts, 
as injectors of chemical elements, as sources of kinetic energy into the interstellar medium and as progenitors of 
supernovae and, likely, as progenitors of long soft 7-ray bursts. 

Evolutionary scenarios, at solar metallicity, for Wolf-rayet stars follow the sequences shown below 
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M > 9OM : 





60- 


90 M : 





40- 


60 M : 





30- 


40 M : 





25- 


30 M Q : 





10- 


25 M : 






- Of - WNL - (WNE) - WCL - WCE - SN (Hypernova low Z ) 

- Of/WNL^LBV - WNL(H poor)- WCL-E - SN(SNIIn?) 

- BSG - LBV WNL -(WNE) - WCL-E - SN(SNIb) 

- WCL-E - WO SN (SNIc) 

- BSG - RSG - LBV ? - WNE - WCE - SN(SNIb) 
-(BSG)- RSG - BSG (blue loop) - RSG - SN(SNIIb, SNIIL) 

- RSG - (Cepheid loop, M < 15 M ) RSG - SN (SNIIL, SNIIp) 



Of-stars are O-type stars showing emission lines. WN stars present emission lines from nitrogen and helium. 
Depending on the line ratios between N III-IV and He I II, they are classified as WNE (WN2-WN5) or WNL (WN7- 
WN9), with WN6 being early (E) or late (L) type. WC spectral type depend on the line ratios of C III and C IV 
(WC4-WC6 are early, WC7-WC9 are late). Oxyg en-rich WO stars extend the WCE sequence, exhibiting strong O VI 
emission lines (see the review b y ICrowtherl . 1 20071 and references therein). LBV stands for Luminous Blue Variables 
(see the review by IVinkl . l2009l and references therein). BSG and RSG for blue and red supergiants respectively. 
The sign «=S> means back and forth motions between the two stages. The various types of supernovae are tentatively 
indicated (see the review bv ISmarttl I2009L for the core-collapse supernova classification). The limits between the 
various scenarios depend on metallicity Z and rotation. The limiting masses indicated corresponds to the standard 
metallicity case, i.e. Z = 0.020. 

From the above evolutionary scenarios, we see that at the standard metallicity, the minimum initial mass for a 
star to enter the WR phase is between 25 and 30 M . This inferior limit is however dependent on the stellar models 
considered and in particular on the mass loss rates used and the mixing processes (size of the convective zone and 
other mixing processes in radiative zones as for instance those induced by rotation). These processes, mass loss and 
mixing, also affect the durations of the WR phases. 

Let us recall some difficulties faced by standard stellar models concerning the WR stars. A good agreement between 
the predictions of the stellar models for the WR/O number ratios and the observed ones at different metallicities in 
regions of constant star formation was achieved provided th e mass loss rates were enhanced by about a factor of two 
during the MS and WNL phases (|Maeder fe Mevnetl . Il994[ ). This solution, which at that time appeared reasonable 



in view of the uncertainties pertaining the mass loss rates, is no longer applicable at present, at least if we take for 
granted the fact that the mass loss r ates during the WR pha se are reduced by a factor 2 to 3, when account is given 
to the clumping effects in the wind (jNugis fe Lamersl |2000( ). Also, the mass lo ss rates for O-type stars have been 
substantially revised (and in general reduced) by the results of IVink et al.l ( 2001). With these revised mass loss rates 



the predicted numbers of WR stars by the standard non-rotating models of Mevnet fc Maederl ( 2003[ ) are much too 
low with respect to the observations. 

A second difficulty of the standard models with mass loss concerns the observed number of transition WN/WC 
stars. These stars show simultaneously some nitrogen characteristic of WN stars and some carbon of the further WC 



stage. The observed frequency of WN/WC stars among WR stars turns around 4.4 % ( van der Huchtj . [200 lh . while 



the frequency predicted b y the standard models without extra-mixing processes are lower by 1-2 orders of magnitude 



( Maeder fc Mevnet] . Il994[ ). This feature cannot be accounted for by a change of the mass loss rates and indicates that 



mixing is not adequately incorporated in the standard stellar models. A third difficulty of the standard models as far 
as WR stars were concerned was tha t there were re latively too many WC stars with respect to WN stars predicted 
at low metallicity (see the review bv lMassevLl2003h . These difficulties are the signs that some process is missing in 
standard models. 



C. Solutions of WR problems in terms of rotation? 

What could be missing? Binarity, changes of the mass loss history, mixing processes can be invoked and there is little 
doubt that playing with these different features, a given observational feature can be reproduced. We think however 
that good stellar models have to reproduce all the main observed features exhibited by massive stars simultaneously 
(as the MS width, the changes of the surface abundances during the MS phase, the variation with the metallicity of the 
of blue to red supergiant ratio, or of the WR/ O number ratio to just cite a few of them) without having to change the 
physics included in order to adapt it to resolve a particular question. This makes the problem much more constrained 
and discard some solutions. As discussed previously, there are some evidences that some extra mixing process is still 
missing in stellar models and that this extra-mixing is probably driven by rotation. Thus in the following we shall 
focus on the consequences for the WR population of rotational mixing. We first indicate some general trends induced 
by rotation and then illustrate the consequences obtained in some specific stellar models. 

Rotational mixing favours the entry into the WR phase in two ways, firstly by allowing chemical species produced 
in the core to diffuse in the radiative envelope and, secondly, by making the mass of the convective core larger. In the 
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non-rotating model, mass loss by stellar winds is the key physical ingredient which allows internal chemical products 
to appear at the surface and thus the formation of a WR star. The star becomes a WR star only when sufficiently 
deep layers are uncovered. In rotating models, the characteristic surface abundances of WR stars can be obtained 
through the effects of mass loss by stellar winds and of rotational mixing. The action of rotation allow WR stars to 
appear through the single star scenario even when the mass loss rates are reduced. To realize that, imagine a star 
rotating so fast that it would follow a homogenous evolution. Such a star can become a WR star, i.e. being a star 
with a log T e ff greater than about 4 and a mass fraction of hydrogen at the surface below ~ 0.4 without losing mass! 

When models with a time-averaged rotation velocity during the MS phase in the range of the observed values are 
considered, then a reasonable number of WR stars can be produced through the single star scenario even using low 
mass loss rates. This is illustrated in Fig. ED Also the rotational models well reproduce the WN/WC ratio at low 
metallicity, the observed fraction of WR stars in the transition stage WN/WC and th e variation with the metallicity 
of the number ratios of type Ibc to type II supernovae (see more detailed discussion in iMevnet fc MaederL I2003L [2005, 
and in Sect. IXI.Aj) . 

The number of WC to WN stars increases with the metallicity (see Fig.[22] 3 ). Man y attempts ha v e been performed 
to reproduce the observed trend: for instance the enhanced mass loss rate models of Mevnet et al.1 ( 19941 ) provided a 
good agreement for solar and higher than solar metallicity but produced too few WN stars in metal-poor regions. The 
inclusion of rotation together with reduced mass loss rates accounting for the effects of clump ing improved the situatio n 
in the metal poor region , but produced too many WN stars at solar and higher metallicities ( Mevnet fc Maederl feoOS'). 
lEldridge k, Vinkl (|2006h show that models that include the mass-loss metallicity scaling during the WR phase closely 
reproduce the observed decrease of the relative population of WC over WN stars at low metallicities. However such 
models severely underestimate the fraction of WR to O-type stars. In that case, to improve the situation, a high 
proportion of Wolf-Rayet stars originating from mass transfer in close binaries would have to be assumed at all 
metallicities. For instance at solar metallicity about 75% of the WR stars should be pro duced in close binary systems 
(see Fig. 5 in lEldridge et all |2008() . This is clearly in contradiction with the results of iFoellmi et al. I (|2003allbh that 
locally in the Milky Way only 24% of the WR stars are binaries. Thus, as shown below, another explanation has to 
be found. 

The WN/WC numbe r ratio depends also on other factors, in particular on the evolutionary scenario. In 
IMevnet fc Maederl (|2003t l2005h . the most massive rotating stars enter into the WR regime already during the MS 
phase. This feature has good and bad effects. On one hand, it allows these models to well reproduce the variation 
of the number fraction of WR to O-type stars since it significantly increases the WR lifetimes. On the other hand, 
it produces very long WN phases since the star enters into the WR phase having still a huge H-rich envelope. As a 
cons equence, too low values for the WC/WN ratio are obtained at solar and higher metallicities. 

In IMevnet fc Maederl (|2003l[2005h . the hypothesis has been made that when a star enters into the WR stage during 
the MS phase, it avoids the Luminous Blue Variable phase. Actually, stars may behave differently. It may well be 
that a star which becomes a WR star during the MS phase, enters a LBV phase after the core H-burning phase, before 
evolving back into the WR regime. When this evolutionary scenario is followed, reasonable values for both the WR/O 
and the WC/WN ratios are obtained. Indeed the ratios of WR/O and of WC/WN given by these models at the solar 
metallicity are 0.06 and 0.9 which compare reasonably well with the observed values of 0.1 and 0.9 respectively. Both 
ratios are not reproduced by the non-rotating models to which a similar scenario is applied. This discussion illustrates 
the possible key role that the LBV phase may play in shaping the WC/WN ratio. 



D. Pop III star models with rotation 

Understanding the evolution of massive stars at low and very low metallicity is a requirement to address questions 
such as the nature of the sources of the reionization in the early Universe, the evolution of the interstellar abundances 
during the early phases of the evolution of galaxies, for finding possible signatures of primordial stellar populations 
in the integrated light of very distant galaxies and for discovering which objects are the progenitors of the long soft 
Gamma Ray Bursts. At present, the most "iron" poor objects known in the Universe are not very far from us since 
th ey are galactic halo field s tars. Provided these stars are trustworthy very metal poor stars (a view recently challenged 
bv lVenn fc Lamberd 1200^ . these objects offer a unique opportunity to study the yields of the first generations of 
stars. 

In this section, we focus the discussion on Pop III stars, i.e. on stars which were made up of material having 



We consider here regions having reached a stationary situation, i.e. regions where the star formation rate can be considered to have 
remained constant for the last twenty million years. 
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been processed only by primordial nucleosynthesis. In an environment with primordial composition, one expects the 
following differences with respect to the more classical evolution at higher metallicity: 

• At strictly Z = 0, the cooling processes, so important for allowing the evacuation of the energy produced when 
the molecular clouds collapse and thus its fragmentation, are not so efficient than when metals are present. 
This favors the formation of mor e massive stars. The in itial mass function is probably different depending on 
the mass of the "minihaloes" (see lGreif fc Bromml l2006h . In minihaloes with masses between 10 6 and 10 s M©, 
virial temperature is between 10 3 and 10 4 K and the cooling is due mainly to the molecules H2. This allows 
the formation of stars with characteristic masses > 100 M©. In minihaloes where ionization occurs prior to 
the late stages of the protostellar accretion process, namely those with a virial temperature superior to 10 4 K 
and thus with masses above 10 s M Q 4 , the hydrogen deuteride (HD) molecule provides an additional cooling 
channel. In those minihaloes, metal-free gaz can cool more efficiently. This leads to the formation of stars with 
masses superior to 10 M©. Thus at the very beginning, we would have first, during a quite short period, only 
very massive Pop III stars and, when gas assembles in more massive haloes (or is reionized by the first stars), 
a second Pop III star generation appears with smaller characteristic masses (about 10 M© and above). When 
the metallicity becomes higher than 10~ 35 Z©, i. e. for Z above about 5 10 -6 a normal IMF governs the mass 
distribution of newly born stars 5 . According to iGreif fc Bromml ( 20061 ) the very massive Pop III stars only 



contribute marginally to feed the reservoir of ionizing photons and to the chemical enrichment of the interstellar 
medium. Much more important are on this respect the less massive Pop III stars born in more massive haloes 
or in relic HII regions. 

• The (nearly) absence of heavy elements implies that massive stars cannot begin burning hydrogen through the 
CNO cycle but through the pp chains. However the energy output extracted from the pp chains is not sufficient 
to compensate for the high luminosity of these stars. The stars must compensate the deficit of nuclear energy 
by extracting energy from the gravitational reservoir, i.e. they contract. At a given point however, due to this 
contraction, the central temperature becomes high enough for activating triple a reactions. Some carbon is then 
produced. When the mass fraction of carbon is of the order of 10 -12 , the CN cycle becomes the main source 
of energy and the evolution then proceeds as in a more metal rich massive stars. For stars above about 20 M©, 
activation of the CN cycle intervenes ve ry early during the c ore H-burning (typically before five percents of the 
hydrogen at the centre is consumed, see lMarigo et al. [ l2003h . 



Pop III stars are more compact due to the contraction they undergo at the beginning of the core H-burning 
phase (see above) and to the fact that the opacity of primordial material is smaller than that of metal rich one. 
Typically a Pop III 20 M© star on the ZAMS has a radius reduced by a factor 3.5 with respect to the radius 
of the corresponding star with Z = 0.020. Even passing from the very low metallicity Z — 10~ 5 to already 
produces a decrease of the radius by a factor 2! Smaller radius favors more efficient mixing in two ways: first, 
for a given value of the diffusion coefficient, D, the mixing timescale decreases with the square of the radius as 
~ R 2 /D, second as can be seen from Fig. [531 m P°P HI models, the meridional currents have smaller amplitudes 
and present a very different configuration than at higher metallicities. Typically at high metallicity (look at the 
curve for Z = 0.020), there is an extended outer cell which tends to reduce the f2-gradient. This is no longer 
the case in metal poor models. 



Recently lEkstrom et al.l ( 2008al) presented a grid of Pop III stellar models including the effects of rotation. Evolu- 



tionary tracks of non-rotating and rotating Pop III stellar models are shown in Fig. [5U (left panel) . The models are 
computed until the end of the core Si-burning, except the 9 M© that has developed a degenerate core before carbon 
ignition and has thus been stopped then, and the 15 M© model that has been stopped at the end of O-burning also 
because of a too degenerate core at that time. An initial velocity of 800 km s _1 on the ZAMS was chosen. For the 
60 M© model, this corresponds to a value of u/w C rit = 0.52 on the ZAMS, which is slightly superior to 0.4, the value 
required at solar metallicity to obtain averaged velocities during the MS phase corresponding to observed values. 

We notice that the ZAMS is shifted toward lower effective temperature and luminosity with respect to the non- 
rotating case as was the case at higher metallicities 6 . Then, when the evolution proceeds, the tracks become more 
luminous, and the main-sequence turn-off is shifted to cooler temperature: the core of the rotating models is refueled 
by fresh H brought by the mixing. It thus grows, leading to an enhancement of the luminosity. 



4 This may also happen in relic HII regions left by the first stars, (sec Johnson & Bromm, 2006). 

5 Some authors argu e that star forma tion switches to more classical mode already when Z = 10~ 6 due to dust production in the early 
Universe, (see Schneider et al., 200j|). 

6 We recall that this shift is due to the sustaining effect of the rotation: the gravity is counter-balanced both by the gas pressure and the 
centrifugal force in such a way that the star behaves like a lower mass one. 
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The onset of the CNO cycle described above can be seen in the HRD: the tracks evolve toward the blue side of the 
diagram, until the energy provided by the CNO cycle stops the contraction and bends the tracks back in the usual 
MS feature. In the rotating 9 M© model, this happens at an age of 12.2 Myr (when the central H mass fraction is 
X c = 0.439) while in the non-rotating one it happens a little earlier, at an age of 10.9 Myr (but at a similar burning 
stage: X c — 0.439). In the case of the non-rotating 15 M© model, it happens after merely 1.5 Myr (X c — 0.695), 
while it takes 2 Myr (X c = 0.677) in the case of the rotating one. 

After central H exhaustion, the core He-burning phase (CHeB) starts right away: the core was already hot enough 
to burn a little He during the MS and does not need to contract much further. This prevents the models to start a 
redward evolution, so they remain in the blue part of the HRD at the beginning of CHeB. Then, something particular 
happens to the rotating models: because of rotational mixing, some carbon produced in the core is diffused toward the 
H-burning shell, allowing a sudden ignition of the CNO-cycle in the shell. This boost of the shell leads to a retraction 
of the convective core and a decrease of the luminosity. At the same time, it transforms the quiet radiative H-burning 
shell into an active convective one. Some primary nitrogen is produced (see more on that point in the next section). 

All the models, except the 9 M Q , reach the critical velocity during the MS phase (see the right panel of Fig. [24]) . 
Once at critical limit, all the models remain at the critical limit until the end of the Main-Sequence phase. On the right 
panel of Fig. [Ml the 85 and especially the 200 M© seem to depart from fi/Q cr it = 1, but this is due to the limit shown 
here being only the f2-limit, where the centrifugal force alone is taken into account to counterbalance the gravity. In 
the two above models however, the radiative acceleration is strong and the mo dels reach the so-ca l led £}r -limit, that 
is the second root of the equation giving the critical velocity: <feg [1 — T] = (jMaeder fc Mevneu l2000bT) . The true 
critical velocity is lowered by the radiative acceleration, and though the f2/fi cr i t ratio plotted becomes lower than 1, 
these models are actually at the ftT critical limit and remain at this limit till the end of the MS phase. At central 
H exhaustion, the moderate inflation of the radius brings the surface back to sub-critical velocities. The mass which 
is lost by the mechanical winds amounts only to a few percents of the initial stellar mass and thus does not much 
affect neither their evolution, nor their nucleosynthetic outputs. Much more mass can be lost by mechanical mass loss 
(see ne xt subsection) when the effects of magnetic fields are accounted for as prescribed in the Tayler-Spruit dynamo 
theory ( Spruitl . 12002 ). or when the metallicity is non-zero. 



E. Strong mass loss in Pop III stars? 



According to iHeger et al.l (120031) . the fate of single stars depends on their He-core mass (M a ) at the end of the 



evolution. They have shown that at very low metallicity, the stars having 64 M© < M a < 133 M Q will undergo pair- 
instability and be entirely disrupted by the subsequent supernova. This mass range in M a has been related to the 
initial mass the star must have on the main sequence (MS) through standard evolution models: 140 M© < M; n ; < 260 
Mq. However the link bet ween the masse s of the He-core mass and the initial mass can be very different depending 
on the physics considered. ^Ek strom et al. | (l2008al) showed that a Z=0, 150 M© stellar model, having an initial ratio 



between the equatorial veloc ity and the c ritical one equal to t>i n i/w C rit = 0-56, c omputed acco unting for the Tayler- 
Spruit dynamo mechanism ( Spruitl . l2002f) and the effects of wind anisotropy ( Maederl . Il999l ) will lose such great 



amount of mass that it will reach the end of the core He-burning phase with a mass of M a too small to go through a 
pair-instability process. 

In Fig. [25l we present the evolution in the HR diagram ( left panel) and the evolution of mass with time (right 
panel) for the 150 M© stellar model of lEkstrom et al. I (|2008aT> . The grey line shows a non-rotating model computed 



with the same physics for comparison. During its whole evolution up to the end of core He-burning, the non-rotating 
model loses only 1.37 M©. This illustrates the weakness of radiative winds at Z = 0. The evolution of the rotating 
model (black line) can be described by four distinct stages: 

1. (continuous part, lower left corner) The model starts its evolution on the MS with only radiative winds, losing 
only a little more than 0.002 M©. During this stage, the ratio of the surface velocity to the critical one increases 
quickly, mainly because of the strong coupling exerted by the magnetic fields. 

2. (dashed part) When the central content of hydrogen is still about 0.58 in mass fraction, the star reaches the 
critical velocity and starts losing mass by mechanical mass loss. It remains at the critical limit through the 
whole MS, but the mechanical wind removes only the most superficial layers that have become unbound, and 
less than 10% of the initial mass is lost at that stage (11.44 M©). The model becomes also extremely luminous, 
and reaches the Eddington limit when 10% of hydrogen remains in the core. Precisely, it is the so-called fir-limit 
that is reached here. 



3. 



( dotted part ) The combustion of helium begins as soon as the hydrogen is exhausted in the core, then the radiative 
H-burning shell undergoes a CNO flash, setting the model on its redward journey. The model remains at the 
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fir-limit and loses a huge amount of mass. The strong magnetic coupling keeps bringing angular momentum to 
the surface and even the heavy mass loss is not able to let the model evolve away from the critical limit. The 
mass lost during that stage amounts to 53.46 M Q . When the model starts a blue hook in the HR diagram, its 
surface conditions become those of a WR star (X sur f < 0.4 and T e s > lO'OOO K). The luminosity drops and 
takes the model away from the r-limit, marking the end of that stage. 

4. (continuous part) The rest of the core He-burning is spent in the WR conditions. The mass loss is strong but 
less than in the previous stage: another 26.34 M Q are lost. 

At the end of core He-burning, the final mass of the model is only Mfi n = 58 Mq , already below the minimum M a 
needed for PISN (M a > 64 Mq). Note that the contraction of the core after helium exhaustion brings the model back 
to critical velocity, so this value for Mg n must be considered as an upper limit. 

This result shows that a fast rotating Pop III 150 M Q may avoid to explode as a PISN. Of course it is by far 
not certain that the conditions required for such a scenario to occur are met in the first stellar generations but it 
underlines the fact that fast rotatio n may drastically change the picture. 

The yields of PISN is discussed bv lHeeer k Wooslevl (l2005h and are dominated by alpha nuclei. The iron-to-oxygen 



ratio show a strong increase with the initial mass. The PISN yields show a strong decrease beyond the iron group, 
have essentially no s-process and r-process element production. These nucleosynthetic signatures of PISN are not 
observed in the most metal poor halo stars. Is this due to the above scenario? To the fact that the signature was very 
quickly erased by the next generations of stars? 7 Or were such high mass stars not formed? These various hypotheses 
cannot be disentangled at the present time, but the observation of more and more metal-deficient stars will probably 
provide elements of response to these questions. 



XI. SUPERNOVAE AND GAMMA RAY BURSTS 

A. The type lb and Ic supernovae 

Type lb supernovae are core-collapse supern ovae whose sp e ctrum shows no hyd r ogen lines. The spectra of type 
Ic show no hydrogen and helium lines (see e.g. iNomoto et all 11994 IWheeler et al. I. Il987h . The progenitors of these 
core collapse supernovae are thus believed to be stars stripped of their original H-rich envelope for type Ib's and also 
of their He-rich envelope for type Ic's. Progenitors are therefore naked stellar cores as e.g. Wolf-Rayet stars. The 
obs erved frequency at solar metallicity of type Ibc supernovae is about 20% the frequency of type II supernovae (see 
e.g. ICappellaro et all I1999D . which represents a significant fraction of all core collapse supernovae. In four cases, the 
typical spectrum of a typ e Ic supernova has been observed associated with a long soft gamma ray burst (GRB) event 
(IWooslev k Blooml . l2006h . indicating a privileged link between type Ic's and the most powerful supernova explosions 



observed in the Universe. 

Fig. [26] shows t he different superno va types expected for various initial masses and metallicities as given from 
rotating models by iGeorgy et al.l (2009). Type Ic SNe cover a much broader range of initial masses than type lb. 



This trend reflects that most stars which, at a given stage, in their evolution have peeled off their H-rich envelope do 
not stop at that stage. Their evolution drives them beyond, up to the stage where most of their He-rich layers have 
been peeled off. 

One immediately sees that the number fraction of type Ibc to type II supernovae should increase with metallicity. 
This is an effect of the stronger winds at higher metallicities, which favor the formation of H-free stars at the 
end of their evolution. Ve ry interestingly, this expected increase of the Ibc/II ratio has been put in evidence by 
iPrantzos k Boissierl (2003). When one compares the the oretical ratios obtained from single stars with the observed 



ones, a good fit is obtained only for the rotating models ( Georgy et all [20091 iMevnet k Maederl . l2005h 



However such an observed trend can also be reproduced invoking close binary-star evolution (see for example 
lEldridge et all 120081 : iPodsiadlowski et all Il992l : IVanbeveren et all 120071) . In this scenario, the hydrogen- rich envelope 
is removed either through a Roche lobe overflow process or during a common envelope phase, producing a WR star. 
Thus we face here the situation where two very different models (single stars with rotation/close binary evolution with 
mass transfer) are both able to give a reasonable fit to the data. Actually both scenarios probably contribute to the 
observed populations of type Ibc supernovae. However, it would be interesting to know their relative importance and 
how their relative importance changes with the metallicity. It might be that both scenarios predict different behaviors 



7 Maybe the metal-poor stars we observe are enriched by more SNe than we actually think, and the later contributions are masking the 
primordial ones. 
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for the way the frequencies of the type lb and type Ic SNe vary as a function of the metallicity. Predictions of single 
star models f o r the variations with the metallicity of the seperate ratios of type Ib/II and of Ic/II are discussed in 
iGeorgv et all ((2009t >. 



B. The progenitors of long soft Gamma Ray Bursts 

One of the most prom ising model for the progenitors of the long soft Gamma Ray Bursts (GRB) is the collapsar 
model bv lWooslevi (119931) . In this model, GRB are the product of the evolution of massive stars which, at the end of 



their evolution, produces a fast rotating black hole. In addition, for the GRB to be visible, the massive star should 
have lost its envelope. This model received a strong support when observations showe d that there is a connection 
between so-called long-soft GRB events and type Ic supernovae ( Wooslev fc Blooml . |2006|) . This kind of supernova can 



be produced by both WC or WO stars. In this way, WC and WO st ars are natura l candi dates to be GRB progenitors. 
However, GRBs are primarily found in metal-poor environments ( Modiaz et all [20081 ). while type Ic's SNe appear 



mainly at high metallicity. Many physical reasons have been in voked to explain why GR Bs seem to occur only at low 
metallicities. Among them are the following (see the review bv lWooslev &: Blooml . 120061) : 



At low metallicity, stellar winds (even during the WR phases) are weaker, thus bringing away small quantities 
of angular momentum. Black holes are also more easily formed since higher final masses are obtained. 

As already emphasized previously, in absence of any magnetic coupling, the transport of the angular momentum 
is mainly due to meridional currents, while the mixing of the elements is mainly driven by shear instabilities. At 
low metallicity the transport of angular momentum between the core and the envelope is less efficient than at 
high metallicities (see Fig. |2"5]) . because of slower meridional currents in metal poor stars. This favors a higher 
angular momentum content of the core when the star is at the end of its nuclear lifetime, it also favors steeper 
f2-gradient in the interior triggering more efficient shear instabilities and thus chemical mixing. 

■ Since the chemical mixing due to rotation is more efficient at low Z, homogeneous evolution is more easily 
obtained in metal-poor regions. Homogeneous evolution allows massive stars to produce a type Ic SN event 
without having to lose large amounts of mass (and thus of angular momentum). Indeed, a perfectly homogeneous 
evolution (actually never realized) would allow the formation of a pure CO core (and then lead to a type Ic SN 
event) at the end of the core He-burning phase without the need for the star to lose any mass! Note that such 
homogeneous evolut ion allows massive stars to produce much more ionizing photons than standard one (see e.g. 
IMevnet et ail . 120081) . 

The distributi on of initial velocities a t low metallicity might contain more fast rotators than at high metallicities 
(see Fig. 9 in iMartavan et al. 1. I2007D . 



It is interesting to compare the observed GRB frequency with the observed frequency of potential candidates. First, 
as mentioned above, type Ic supernovae at low metallicity do appear interesting candidates. From Fig. 1271 one can see 
that the observed rate of type Ic supernovae from single star models is still above the estimated number ratio GRB 
/ core collapse supernovae (CCSNe) even when one only considers low-metallicity, type Ic SNe. Here it is supposed 
that the formation of a BH does not prevent a SN event. In case, only rare circumstances would allow a SN to occur 
when a black hole is formed, then the situation might be very different. 

Other interesting candidates are the WO stars that primarily occur at low metallicity. These stars will explode as 
a type Ic SN. The frequency of type Ic's SNe with WO star progenitors is shown in Fig. [27] The expected rate is only 
marginally compatible with the observed GRB rate (assumi ng that the aperture angle of the bipolar jet is very small, 
typically around 1°). This conclusion has been obtained bv lHirschi et al.l (2005a), so even restraining the progenitors 



of GRB to WO stars would still not match the observed frequency of GRBs. 

From the above discussion, it can be deduced that exploding as a type Ic SN in metal-poor regions (or having a 
WO progenitor) is not a sufficient condition for obtaining a GRB. Physical characteristics shared by a subsample of 
the metal poor type Ic events exists that are needed to obtain a GRB event. Proba bly this physical characteristic is 
the high angular momentum in the core (IWooslev fc Heeerl . 120061: lYoon et all 120061 ) . 



If GRB would only occur for initial l y very fast-rotating star s , rotating so fast that these stars would follow a 
homogeneous evolution (|Maederl . Il987al: iMevnet fc Maedert 120071: IWooslev & Hegerl 120061: lYoon et all 1 20061) . can we 



expect to find any peculiar feature in the chemical composition of the ejecta testifying this previous homogeneous 
evolution? Or in other words, is there any difference in the chemical composition of the ejecta between a type Ic 
having "normal WC or WO" progenitors and those arising from a mode l that followed a ho mogeneous evolution 
during the MS phase? The answer is probably no. It has been shown by IGeorgv et al. (2009) that the differences 
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in the masses of He, CNO elements, and Z are very small. Thus there is little chance from the observations of the 
composition of the ejecta in CNO elements to be able to distinguish between a normal and a homogeneous evolution. 

A point that would be interesting to check is the following: since the star rotates very fast during a great part of 
its evolution, it will produce anisotropic stell ar winds (see F ig. |2"5]) . typically bipolar winds and probably equatorial 
mass loss when the crit ical limit is reached (jMaederl . 19991) . These anisotropic winds will shape the circumstellar 
environment of the star ( van Marie et al 1. 120081 and see also Fig. [28]) and it might be that some traces of the resulting 
particular morphology will still be present at the time of the SN event (for instance, equatorial mass loss probably 
gives rise to a slow equatorial expanding disk whose traces might still be present when the star explodes). In that 
case the circumstellar environment of GRBs may be peculiar. This can in turn have an impact on some features in 
the spectrum. More generally, the circumstellar environment of stars that have been in a not too remote past very 
fast rotators may also be characterized by such features. 



XII. STELLAR NUCLEOSYNTHESIS 
A. Rotation and stellar yields at Z=0.020 

Figure [23] displays the total stellar yields (i.e. the mass of an element newly synthesized by the star and ejected 
into the interstellar medium) divided by the init ial mass of the star, p *^, as a function of its initial mass, m, for the 
non-rotating (left) and rotating (right) models ( Hirschi et all l2005b[ ). The different shaded areas correspond from 



top to bottom to for 4 He, 12 C, 16 and the rest of the heavy elements. The fraction of the stars locked in the 
remnant as well as the expected explosion type are shown at the bottom. The dotted areas show the wind contribution 
for 4 He, 12 C and 16 0. 

The remnant masses have been estimated using the rela tion between the mass of th e carbon-oxygen core and the 
mass of the remnant given bv lMaederl ( 19921 ) (see details in lHirschi et al. I. l2004l2005bh . This value is quite uncertain 



and may affect the yields, especially those of the iron-group elements and heavier. The yields shown in Fig. [25] arc 
those obtained before explosive nucleosynthesis occurs. 

For 4 He (and other H-burning products like 14 N), the wind contribution increases with mass and dominates for 
M > 22 M Q for rotating stars and M > 35 M for non-rotating stars, i. e. for the stars which enter the WR stage. 
For very massive stars, the SN contribution is negative and this is why P4^ cra is smaller than p™^- In order to 
eject He-burning products, a star must not only become a WR star but must also become a WC star. Therefore for 
12 C, the wind contributions only start to be significant above the following approximative mass limits: 30 and 45 M 
for rotating and non-rotating models respectively. Above these mass limits, the contribution from the wind and the 
pre-SN are of similar importance. Since at solar metallicity, no WO star is produced, for 16 0, as for heavier elements, 
the wind contribution remains very small in the whole mass range. 

So we see that, below ~ 30 M , rotation increases the total metal yields, Z, and in particular the yields of carbon 
and oxygen by a factor of 1.5-2.5. As a rule of thumb, the yields of a rotating 20 M star are similar to the yields of 
a non-rotating 30 M star, at least for the light elements shown in Fig. [2S] We see therefore that at Z — 0.020, the 
effects of rotation are qualitatively similar to those obtained by models with larger convective cores as for instance 
those computed with an overshoot. This is of course not surprising since rotational mixing tends also to increase the 
core size at a given evolutionary phase. For very massive stars (~ 60 M ), rotation increases the yield of helium but 
does not significantly affect the yields of heavy elements. The similarities of the results between the rotating and 
non-rotating models for the 60 M model comes from the fact that, although the mass is not lost at the same stages 
in both series of models, the total mass loss is nevertheless not very different. The non-rotating model goes through 
an LBV phase after the MS phase where a significant fraction of the mass is lost. The rotating model enters into the 
WR phase already during the MS phase and no LBV phase has been accounted. 

Figure |3"01 presents the stellar yields convolved with the Salpeter initial mass function (IMF) (dN/dM oc M~ 2 - 35 ). 
This reduces the importance of the very massive stars. Nevertheless, the differences between rotating and non-rotating 
models remain significant, especially around 20 M . Despite these differences, it will be difficult to discriminate 
between these two series of yields from observations of the evolution of the abundances as a function of the increasing 
metallicity, because in such comparisons, the yields are only one of the ingredients of the chemical evolution models. 
Other factors as the star formation history, the initial mass function, the rate of infall and outflow etc... will blurr the 
picture. More constraining are the early phases of the chemical evolution of the galaxies during which only massive 
stars had time to contribute. 



30 



B. Impact on the early chemical evolution of galaxies 

Very metal poor halo stars hav e formed (at least in part) from matter ejected by very metal poor massive stars (see 
e.g. [Chiappini et all . I2008L 12006ft . By very metal poor stars we mean here Pop III stars and/or stars with an initial 
metal abundance lower or equal to that of the metal poor halo considered. Their surface compositions thus reflect 
the nucleosynthesis occurring in the first generations of massive stars (provided of course that no other processes as 
accretion or in-situ mixing mechanism has changed their surface composition) . 

Interestingly, many observations of these stars show puzzling features . Among them let us cite the two following 
ones shown in Fig. |32] spectroscopic observations (e.g. ISpite et al 1 12005ft indicate a primary production of nitrogen 
over a large metallicity range; Halo stars with log(Q/H)+12 inferior to about 6.5 prese nt higher C/O ratios than halo 
stars with log(0/H)+12 between 6.5 and 8.2 (jAkerman et al.1 . 12004 ISpite et al.1 . |2005( ). 



Fast rotating massive stars are very interesting candidates for producing primary nitrogen at low metallicity and 
in a very short time delay as is required by the observations of the N/O plateau shown in the upper panel of Fig. 1321 
Primary nitrogen production o ccurs in both high mass an d intermediate mass stars provided they rotate and have a 
metallicity below about 0.001 (jMevnet fe Maederl . [2b02allbh . The initial rotation needed to obtain sufficient primary 
nitrogen production corresponds to values between 50 and 70% of the critical rotation on the ZAMS. These values are 
above the usual ratio adopted for solar metallicity models, which is 40%, but are still far from extreme values near 
the critical limit. Let us note that recent models of Pop III star formation seem to support very high initial rotational 
velocities (|Stacv et al.L l2010h . 

For a given initial velocity, primary nitrogen production is more efficient when a small amount of metals is present 
than in stars made of primordial material. Indeed, although rotational mixing in Pop III stars is by far not a negligible 
effect, it remains at a relatively modest level due to the absence of strong contraction at the end of the core H-burning 
phase. On the contrary, when Z >~ 10~ 10 , the physical conditions during the core H-burning phase and the core 
He-burning phase are so different that a strong contraction occurs at the end of the core H-burning phase leading 
to steep f2- gradient, strong mixing and important primary nitrogen production. For metallicities higher than about 
0.001, rotational mixing is not efficient enough for triggering important primary nitrogen production (at least for the 
rotational velocities corresponding to the observed ones at this metallicity) and thus rotational mixing, although still 
important for explaining the surface enrichments, does not change the stellar yields as much as at very low metallicity. 

A comparison of various yields for CNO elements for Pop III stars and very metal poor stars is shown in Fig. 1311 
For 12 C (left) and 16 Q ( right), our non-rotating 15 and 25 M Q yields compare within a factor 2 or 3 with those of 



IChieffi fc Limongil ( 20041 ). Their 25 produces also a huge amount of 14 N (centre), even a factor ~ 10 higher than 



ours. Recent non- rotating models bv lHeger fc Wooslevl ( 2010T ) predict also important quantities of primary nitrogen 



in some mass range, indicating that primary nitrog en production in Z=0 star models may not be linked to rotation. 
However, according to chemical evolution model of IChiappini et all (|2006ft . the production of primary nitrogen only 
by Pop III stars may not be sufficient to account for the observed N/O plateau shown by halo stars. Primary nitrogen 
production over a metallicity range is needed. At the moment, rotating models can provide a mechanism which 
operates over a sufficient extended metallicity range to account for the observations. 

When the yields of rotating massive star models are included in chemical evolution model for the halo, nice 
agreement can be obtained with the observations. It is possible to account for the high observed N/O ratio at 
very low log(0/H)+12 values and to reproduce the observed C/O upturn mentioned just above. This is illustrated in 
Fig.[32l where predictions for the evolution of N/O and C/O of chemical evolution models using different sets of yields 
are compared (Chiappini et al. 2006a 8 ). We see that the observed N/O ratio is much hig her than what is p r edicted 



by a chemical evolution model using the yields of the slow-rotating Z = 10 5 models from Mevnet fc Maeder d20 02a) 



down to Z = 0. When adding the yields of the fast-rotating Z = 10~ 8 models from iHirschil (|2007ft : iMevnet et al.l 
(2006) 9 the fit is much improved. The same improvement is found for the C/O ratio, which presents an upturn at 
low metallicity. These comparisons support fast rotating massive stars as the sources of primary nitrogen in the most 
metal poor galactic halo stars. 



Hig h N/O and the C/O upturn of the low-metallicity stars are also observed in low-metallicity DLAs (|Pettini et al.1 . 
2008, see the crosses in Fig. [32]) . We note that the observed points are below the points for the halo stars in the N/O 
versus O/H plane. This may be attributed to two causes: either the observed N/O ratios observed in halo stars are 
somewhat overestimated or the difference is real and might be due to different star formation histories in the halo 
and in DLAs. Let us just discuss these two possibilities. 



8 The details of the chemical evolution models can be found in Chiappini et al. (2006b), where they show that such a model reproduces 
nicely the metallicity distribution of the Galactic halo . This means that the time scale for the enrichment of the medium is well fitted. 

9 These models use the same physics as the models of Meynct & Macdcr (2002a), simply the initial rotation velocities considered are 
higher. 
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Measures of nitrogen abundances at the surface of very metal poor stars is quite challenging, much more than 
the measure of nitrogen in the interstellar medium as is done for the DLAs, therefore one expects that the data for 
DLAs suffer much smaller uncertainties than those for halo stars. In that respect the observed N/O ratios in DLAs 
should be preferred to compare with the chemical evolution models. In case both DLAs and the galactic halo had 
the same star formation history, i.e. are the result of an intense and rapid star formation epsiode, then the chemical 
evolution models presented in Fig. [32] can also apply to DLAs. We see that the DLA data still require models with 
fast rotation to be fitted (the shift of the DLA data towards lower N/O ratios is by far not sufficient to discard fast 
rotating models). 

Probably the star formation history in DLAs is not the same as in the halo. While in the halo we see the result of 
strong and rapid star formation episodes, in DLAs one might see the result of much slower and weaker star formation 
episodes. In that case, both massive stars and intermediate mass stars contributed to the build up of the chemical 
abundances and the chemical evolution models presented in Fig. [32] do no long apply to theses systems. It will 
be very interesting to study the results of chemical evolution models adapted to this situation and accounting for 
stellar yields from both massive and intermediate mass stars. Let us just mention at this stage that primary nitrogen 
production in meta l poor i ntermediate mass stars is also strongly favoured when rotational mixing is accounted for 
(|Mevnet fc Maederl . l2002aHbh 10 . Thus also in that case, rotation may play a key role. 

The p rimary nitrogen produc tion is accompanied by other interesting features such as the production of primary 
13 C (see IChiappini et al. , 2008). and of primary Ne. Primary Ne is produced by diffusion of primary nitrogen 
from the H-burning shell into the core He-burning zone, or by the engulfment of part of the H-burning shell by the 
growin g He-buning core. These processes occur in rotating massive star models ([Hirschil . l2007t iMevnet fc Maedeil 
l2002aD . In the He-burning zone, 14 N is transformed into 22 Ne through the classical reaction chain 14 N(a,7) 18 F(/3 + 
i/) 18 0(a, 7 ) 22 Ne. 

In the He-burning zones (either in the core at the end of the core He-burning phase or in the He-burning shell during 
the core C-burning phase), neutrons are released through the reaction 22 Ne(a,n) 25 Mg. These neutrons then can either 
be captured by iron seeds (or elements with atomic numbers near the one of the iron called iron peak elements) and 
produce s-process elements or be captured by neutron poisons and thus removed from the flux of neutrons which is 
useful for s-process element nucleosynthesis. From what precedes one can easily conceive that the final outputs of 
s-process elements will depend on at least three factors: the amounts of 1.- 22 Ne, 2.- neutron poisons and 3.- iron 
seeds. In standard models, when the metallicity decreases, the amount of 22 Ne decreases (less neutrons produced), 
the amount of neutron poison becomes relatively more important because for [Fe/H] 11 < -2 main neutron poisons 
are primary elements and the amount of iron seeds decreases also. Thus very small quantities of s-process elements 
are expected (see the triangles in Fig. l33l) . When primary nitrogen and therefore primary 22 Ne is present as given 
by rotating models which can reproduce the observed trends for the N/O and C/O ratios observed in the halo stars, 
then a very different output is obtained. One can see that globaly the abundances of the s-process elements are 
increased by many order of magnitudes, and that the elements produced in the greatest quantities are no longer in 
the iron-nickel region as is the case in standard models but cover an extended region from strontium up to baryum 
(|Pignatari et all I2OO8I) . 

These first results need to be extended for other masses, rotation and metallicities. However they already show 
that some heavy s-process elements, not produced in standard models, might be produced in significant quantities in 
metal poor rotating stellar models. It will be very interesting in the future to find some non ambiguous signature of 
the occurrence of this process in the abundance pattern of very metal poor halo stars. 



C. The Carbon-Enhanced- Metal- Poor stars 

At least 1 star out of 4, with [Fe/H] < -2.5, exhibits [C/Fe] > +1.0 (see the review bv lBeers & Christliebl . l2005h . 
These stars are known as C-enhanced metal poor stars (CEMP). The two most iron poor stars presently known 
([Fe/H] equal to -5.96 and -5.4) are C-rich, they are called C-rich Ultra-Metal Poor Stars (CRUMPS). The origin 
of their peculiar surface abundances is not known. The high nitrogen abundances observed in some CEMP stars 
implies that the material which is responsible for their peculiar abundance pattern must be heavily loaded in primary 
nitrogen. As seen above, rotating stars (both massive and intermediate mass stars) can produce large amounts of 



10 In that paper, the evolution was stopped in the E-AGB phase and was not pursued in the TP-AGB stars, however the quantity of 
primary nitrogen obtained at that stage will be a lower limit of what will be found at the end of the TP-AGB phase and thus represents 
some lower limit value. 

11 [Fe/H]=lg(n(_Fe)»/n(H)») — lg(n(Fe)Q /n(H)Q), where n(Fe)* is the number density of iron nuclei in the considered star and n(Fe)Q 
in the Sun, n(H) t and u(H)q the number of hydrogen nuclei in the star and in the Sun. 
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primary nitro g en. Rotatin g models can account also for the observed enhancements in C, O, Na, Mg and Al, thus 
IMevnet et al.l ( 2006T ) and iHirschil (2007) have proposed that CEMP stars may be formed from material ejected by 
fast rotating stars (massive or intermediate mass stars). Different chemical signatures are expected depending on the 
fact that the CEMP stars are formed from wind ejecta, from wind ejecta together with supernova ejecta or from the 
envelope of an AGB star. For instance CEMP stars formed from wind material of massive rotating stars mixed with 
small amount of pris tine interstellar medium produce stars which are He-rich, Li-depleted and with low 12 C/ 13 C ratios 
IMevnet et al.l ( 2010h . By He-rich we mean stars that, at very low [Fe/H], present a helium mass fraction between 0.30 
and 0.60. 

Rotation is probably a key ingredient to explain the abundance patterns of CEMP stars. Similar non-rotating 
models, without any other extra-mixing, do not succeed to simultaneously explain the enhancements in the three 
CNO elements. 



D. The chemical anomalies in globular clusters 

The finding of a double sequence in the globular cluster u> Cen, as well as in many other globulars, and the further 
interpreta tion of t he blu er sequence by a strong excess of helium constitutes a major enigma for stellar and galactic 
evolution (|NorrisLl200l . The bluer sequence with a metallicity [Fe/ H]= -1.2 or Z = 2-10 3 implies an He-content 
Y = 0.38 (0.35-0.45), i.e. an helium enrichment AY = 0. 14 (see also lNorrisL l2004h. In turn, this demands a relative 



helium to metal enrichment AY/AZ of the order of 70 (jPiotto et al 



2005 ). This is e normous and more than one 



order of magnitude larger than the current value of AY/AZ = 4 — 5 ( Pagel et al.lll992l) obtaine d from extragalactic 
HII regions. A value of 4 -5 is consistent with the ch emical yields from supernovae (|Maederl . ll992l ) forming black holes 
above about 20-25 M Q . iMaeder &: Mevnetl ( 20061 ) showed that the wind contributions of low Z rotating stars are 



able to produce the high AY/AZ observed in the blue Main Sequence in uj Cen. 

Recent observations show that He-r ich stars in globula r clust ers might be a common feature. For instance, clusters 
like NGC 2808, M13 and NGC 6441 (jCaloi fc D'Antonal . 120071 ) show a well-developed blue horizontal branch, and a 



strong slope upward from the red clump to the blue of t he RR Lyrae regions. Both features could be explained if 
He-rich stars are present in these clusters. iKavirai et al.l (|2007l) studied the UV and optical properties of 38 massive 
globular clusters in M87 (elliptical galaxy in virgo cluster). A majority of these clusters appear extremely bright 
in the far-UV. Canonical models (i.e. models without any super He-rich stars) would imply ages for these clusters 
about 3-5 Gyr larger than the age of the Universe according to the WMAP data! This difficulty is removed when a 
super- He-rich (AY/AZ >~ 90) stellar component is supposed to be present in these clusters. 

It has also long been known that globular cluster stars present some striking anomalies in their content in other 
light elements 12 : while in all the Galactic globular clusters studied so far one finds "normal" stars with detailed 
chemical composition similar to that of field stars of same metallicity (i.e., same [Fe/H]), one also observes numerous 
"anomalous" main sequence and red giant stars that are simultaneously d eficient (to various degrees) in C, O, and 
Mg, and enriched in N, Na, and Al (see the review bv lGratton et al. 1. 12004) . 

It is clear now that these chemical peculiarities have been inherited at birth by the low-mass stars we observe today 
in globular clusters. The observed abundance pattern indicates that the material from which these stars formed was 
nuclearly processed by CNO burning and the NeNa- and MgAl-chains that occurred in an early generation of more 
massive and faster evolving globular cluster stars. Therefore, there were at least two generations of stars in Galactic 
globular clusters. The first one corresponds to the bulk of "normal" stars born with the pristine composition of the 
protocluster gas; these objects are those with the highest O and Mg and the lowest He and Na and Al abundances 
also found in their field contemporaries. The second generation contains the stars born out of material polluted to 
various degrees by the ejecta of more massive stars, and which present lower O and Mg and higher He, Na and Al 
abundances than their first generation counterparts. 

Any model aiming at explaining the chemical properties of globular cluster stars should at least give an answer to 
the following questions: (1) Which type of stars did produce the material enriched in H-burning products? (2) What is 
the physical mechanism responsible for selecting only material bearing the signatures of H-processing? (3) Why does 
this process occur onl y in globular clusters? 13 Mod els for explaining t he chemical anomalie s in globular clusters invoke 
either TP- AGB stars (I Ventura fc D'Antonal . 120081) , massive binaries (|de Mink et al.l . l2009h and massive rotating stars 



( Decressin et al.l . l2007bl) . These last authors propose that the matter from which the stars rich in H-burning products 



On the contrary, the content in heavy elements is fairly constant from star to star in any well-studied individual Galactic globular cluster 
(with the notable exception of ui Cen). 

Indeed up to now the peculiar chemical patterns observed in globular clusters, i.e., the O-Na and Mg-Al anticorrelations, have not been 
found in field stars. 



33 



are formed, has been released by slow winds of fast rotating massive stars. Of course, part of the needed material can 
also be released by AGB stars. The massive star origin presents however some advantages: first a massive star can 
induce star formation in its surrounding, thus two effects, the enrichment and the star formation can be triggered by 
the same cause. Second, t he ma ssive star scenario allows to use a less flat IMF than the scenario invoking AGB stars 
(jPrantzos fc Charbonnel 120061) . The slope of the IMF might be even a Salp e ter's o ne in case the globular cluster 
lost a great part of its first generation stars by tidal stripping ( Decressin et all l2007al ) . Finally the material released 
through the mechanism invoked above presents the property that the sum of the C+N+O elements remain constant. 
Interestingly, the sum [(C+N+0)/Fe] for individual globular cluster stars presenting otherwise very different chemical 
features is consta nt to within exp e rimen tal errors in all the clusters studied so far with the exception of NGC 1851 
(see references in iDecressin et all . [2009). These last authors show also that rotating AGB models lead to a large 
increase of the total C+N+O value, which disfavors these models for explaning the origin of the chemical anomalies 
in most of the globular clusters. 



XIII. CONCLUSIONS 

The evolution of massive stars depends sensitively on the metallicity which has an impact on the intensity of the 
line driven stellar winds and on rotational mixing. We can distinguish four metallicity regimes: 1.- the very metal 
poor regime < Z <^ 10~ 10 ; 2.- The low metallicity regime 10 -10 < Z < 0.001; 3.- The near solar metallicity regime 
0.001 < Z < 0.020; 4.- The high metallicity regime 0.020 < Z. In each of these metallicity ranges, some specific 
physical processes occur. 

In the (likely rare) Pop III regime, lEkstrom et al. (2008a) find that models rotate with an internal profile f2(r) close 



to local angular momentum conservation, because of a very weak core-envelope coupling. Rotational mixing drives 
a H-shell boost due to a sudden onset of CNO cycle in the shell. This boost leads to a high 14 N production, which 
can be as much as 10 6 times higher than the production of the non-rotating models. Due to the low mass loss and 
the weak coupling, the core retains a high angular momentum at the end of the evolution. The high rotation rate at 
death probably leads to a much stronger explosion than previously expected, changing the fate of the models. 

In the low metallicity regime, rotating models produce large amounts of primary 13 C, 14 N and 22 Ne. Rotation has 
also probably a strong impact on the s-process elements nucleosynthesis. Rotational mixing, by increasing the CNO 
surface abundances, might trigger important mass losses through radiatively driven stellar winds. Stars may also lose 
some mass by mechanical mass loss when rotating at the critical limit. For all these reasons, rotation may play an 
important role in explaining the CEMP stars and the chemical anomalies in globular clusters. 

In the near solar metallicity regime, rotation and mass loss by stellar winds are of similar importance, however 
the balance of the two effects varies with mass. Neither of the two aspects can be neglected. This metallicity range 
is also the one in which models can be checked and calibrated by comparisons with well observed features either of 
individual stars or of stellar populations. Among these observed features let us cite 

• The observed changes of the surface abundances. 

• The observed changes of the surface velocities. 

• The width of the Main Sequence band. 

• The shape of fast rotating stars measured by interferometric technics, the shape of the nebulosities around high 
mass losing stars, the anisotropies of the stellar winds. 

• In a near future, the size of the convective core as deduced from asteroseismic analysis, as well as the varia- 
tion with depth of the angular velocity as obtained from asteroseismic observations. Note that simultaneous 



observations in asteroseismology an d interferometry (jCunha et all 120071 ). together with estimate of the surface 



magnetic field ( Wade et all 1201 Oal) will provide new and important clues on massive star evolution 



• The existence and variation with Z of the populations of Be stars. 

• The variation with Z of the blue to red supergiant ratios. 

• The variation with Z of the Wolf-Rayet populations. 

• The rotation rates of young pulsars. 

• The variation with Z of different core collapse supernova types. 
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As a general statement, it does appear that models including the effects of rotation provide a much better fit to most 
of the above observed features. Other aspects, not discussed in the present review as the synthesis of 26 Al, the origin 
of the isotopic anomalies in the galactic cosmi c rays and in meteorites can also be affected by rotation ( Arnould et all 
l2006t iBinns et all 120051: iPalacios et ail . 120051) . 

Above solar metallicities, radiative line driven winds become the dominant factor affecting the evolution of massive 
stars. In models with moderate rotation, we note however that the effects of rotational mixing are still important 
but their impact is less apparent being somewhat mixed with those of the stellar winds. An interesting difference 
which occurs at high metallicity is the evolution of the angular momentum. One expects that at high metallicity, 
angular momentum is more easily transported from the core to the surface, due to faster meridional currents, and 
more easily ejected at the surface by the stellar winds. Thus, everything being equal, one would expect that the 
angular momentum of the central regions will be lower when the metallicity is higher. Together with the fact that 
Black-Hole formation is probably more difficult at high Z (also because of the strong stellar winds), this makes the 
formation of collapsars which are considered as a serious candidate for long soft gamma ray burst much less favorable 
at high metallicity. In case the rotation rate of young pulsar depends to some extent to the rotation rate of the core 
at the presupernova stage, the above line of reasoning would lead to the conclusion that the rotation rate of young 
pulsars should be slower in metal rich regions than in metal poor ones. 

As we can see, rotation and mass loss are two key physical processes which have a strong impact on the evolution 
of the massive stars during the whole cosmic history. These two processes, which interact with each other, are 
responsible for a great diversity of outputs, going from the shape of the nebulae, to the evolution of the abundances 
in the interstellar medium, passing through the origin of the Gamma Ray Bursts. It is interesting to note that the 
energy in form of rotational energy in stars is quite small, at most of the order of the percent of the gravitational 
energy. Similarly, metallicity or heavy elements, which shape the intensity of the line driven winds, represent only at 
most two percents of the mass of all chemical elements. Nevertheless, the impacts of both rotation and metallicity 
are very important. It reminds us that sometimes a pinch is sufficient for creating new worlds! 
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FIG. 1 Pro bability density by km s of rotation velocities for 496 stars with types 09.5 to B8, i.e. masses between about 3 
and 20 M (|Huang fc Giesl , l2006al ). 
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FIG. 2 Shape of the surface for various values of u 



S2 crit 



(labelled at the bottom of each curve). The x-axis is the equatorial 



radiu s, and the y-axis the polar one. Hence this is how we would see the star equator-on. R p is the polar radius. Figure taken 
from lGeorgy et all i|201ll ). 
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FIG. 4 The mass fluxes around a rotating star of 100 Mq with 10 6 5 and a ratio of the angular velocity to the break-up 
angular velocity cj = 0.80, assuming a polar T e a = 30000 K ( Mae der" fc Desiacquesl . l200ll ) . 
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FIG. 5 Schematic structure with stream lines of meridional circulation in a rotating 20 Mq model of 5.2 R.0 with Y — 0.30, 
Z — 0.02 and Vim = 300 km s -1 at the beginning of the MS phase. The figure is made as a function of M r . In the upper 
hemisphere on the right section, matter is turning counterclockwise alo ng the outer stream line an d clockwise along the inner 
one. The inner sphere is the convective core. It has a radius of 1.7 R© (Meynet & Macdcr, 2002a). 
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FIG. 6 Unstable displacements with m — 1 of the azimuth al magnetic fi eld in the polar region. The large arcs (thin broken 
lines) indicate the horizontal stellar surfaces. Adapted from lSpruitl (|l999f ). 
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FIG. 7 Left: evolution of the angular velocity as a function of the distance to the center in a 20 Mq star with vi u i = 300 
km s" 1 . X c is the hydrogen mass fraction at the cente r. The dotted line shows the profile when the He-core contracts at the 
end of the H-burning phase ijMevnet &i Maederl . l200(jh . Right: rotation profiles at various stages of evolution (labeled by the 
central H content X c ) of a 15 Mq model with X = .705, Z = 0.02, an initial velocity of 300 km s _1 and magnetic field from 
the Tayler-Spruit dynamo l|Maeder fc Mevnetl l2005bT ). 
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t(10 6 y) 



FIG. 8 Evolutio n with mass loss (only) of the total mass and internal structure as a function of time for a star with an initial 
mass of 60 M l|Maederj. Il987af) . The cloudy regions indicate convection, the heavy diagonal hatching the regions of high 
nuclear production, the light vertical hatching regions of variable composition. 




FIG. 9 Variation of the diffusion coefficients as a function of the radius inside a 20 Mq model with an initial metallicity Z=0.014 
and an initial velocity on the ZAMS equal to 40% the critical velocity. From left to right, is shown the situation corresponding 
to central mass fraction of hydrogen equal to 0.60, 0.30 and 0.01 respectively. Courtesy of S. Ekstrom. 
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FIG. 10 The differences in log(N/H) as a function of the initial masses at 3 stages during the MS phase for models at Z = 0.02 
with the average rotation velocities (i.e. 217, 209, 197, 183, 172, 16 8 km s" 1 for respect ively 12, 15, 20, 25 40 and 60 M ). 
The 3 stages are indicated by the value of the central H-content X c (|Maeder et al.l . l2008l ). 
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FIG. 11 The evolution of the differences in log(N/H) during the MS phase as a function of the actual rotatio n velocities (sini 
being equal to 1 here) for models of 20 Mq with Z = 0.02 and different initial velocities (|Maeder et al.l . l2008l ). 
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FIG. 12 Evolutionary tracks in the plane surface N/C ratio, normalized to its initial value, versus the rotational period in 
hours for different initial mass stars, various initial velocities and for the metallicities Z—0.02 and 0.002. Positions of some 
periods in days are indicated at the bottom of the figure. The dotted tracks never reach the critical limit during the MS phase. 
The short dashed tracks reach the critical limit during the MS phase. The dividing line between the shaded and non-shaded 
areas corresponds to the entrance into the stage when the star is at the critical limit during the MS evolution. Big circles along 
some tracks indicate the stage when v/v CI i t becomes superior to 0.7. If Be stars are rotating at velocities superior to 70% of 
the critical velocity, present models would predict that they would lie in the region comprised between the big circles and the 
dividi ng line. Beware the d ifferent vertical scales used when comparing similar masses at different metallicities. Figure taken 
from lEkstrom et all (|2008bT ). 
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v sin i (km s _1 ) 

FIG. 13 The N abundance (in a scale whe re log H = 12 . 0) as a function of usini for the MS stars (black dots) in Nil with 
masses between 14 and 20 Mq according to Hu nter et al.l (|2008| ). The downward arrows indicate abundance upper limits. The 
binaries are shown by a square. The evolved stars in a band of 0.1 dex in \ogT e g beyo nd the end of the MS are shown with 
open symbols. The grey band indicates uncertainties of ±0.25 dex (Macd er et alj . l2008h . 
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FIG. 16 Angular velocity a s a function of the l agrangian mass coordinate, m r inside the 25 Mq model (t?i n i= 300 kms 1 ) at 
various evolutionary stages (jHirschi et all . 120041 ). 
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FIG. 18 Evolutionary tracks for non-rotating (dotted lines) and rotating (continuous lines) models with solar metallicity. The 
rotating models have an initial velocity «i n i of 300 km s - . For purpose of clarity, only the first part of the tracks for the 
most massive stars (M > 40 Mq) is shown. Portions of the evolution during the WR phase for the rotating massive stars are 
indicated by short-dashed lines. The long-dashed track for the 60 Mq model corresponds to a very fast rotating star (i>i n i 
~ 400 km s -1 ), wh i ch fol lows a nearly homogeneous evolution. Only the beginning of its evolution is shown. Figure taken from 
iMevnet fc Maedeij (|2000H . 
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FIG. 20 Comparison of the internal distribution of helium in two models of 20 Mq at the middle of the He-burning phase. 
The dashed-dot line c oncerns the models with ze ro rotation and the continuous line represents the case with Uj n j = 300 km 
s~ . Figure taken from iMaeder fc Mevnetl l|2001bh . 
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FIG. 21 Variation of the number ratios of Wq lf-Rayet stars to O-type stars as a function of the meta llicity. The observed 
points are taken as in lMaeder fc Mevnetl (| 19941 ). The dotted line shows the predictions of the models of iMevnet et al. I (|!994h 
with normal mass loss rates. The continuous and the d ashed lines show the predictions of the rotating and non-rotating stellar 
models respectively of IMevnet fc Macdcr (2003, 2005|). The black pentagon shows the ratio pre dicted by Z=0.040 mo dels 
computed with the metallicity dependence of the mass loss rates during the WR phase proposed bv lCrowther et all l|2002h . 
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FIG. 22 Variation of the number ratios of WN to WC stars as a function of metallicity. The grey area encompasses the 
observed ratios. Individ ual measures are indicated by blac k circles labeled with the name of the Galaxy (see references in 
IMevnet fc Maedd. 12005ft. So l ar (O /H) value is taken from lAsplund et all l|2005ft , the (O/H) values for the SMC and LMC 
are taken from 111 1 mil <-r el all, 120071 ). The dotted lines show the predictions of the rotating and non-rotating stellar models 
of IMevnet fc Maederl (|2003l . |2005| ). The black pentagon shows the ratio predicted by Z=0.040 models computed with the 
metallicity dependence of the mass loss rates during the WR phase. The open triangle s hows the WC/WN ratio obtained from 
the rotating models allowing an LBV phase during the WR phase ijMevnet et al l l200Sft . 
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FIG. 23 Internal profile of U(r) in models of 20 Mq at various metallicities, with f2int/f2 C rit = 0.5, where u(r,8) the vertical 
component of the velocity of the meridional circulation is u(r, 9) = U(r) P2(cos6 l ) with P2 the second Legendre polynomial. 
The radius is normalized to the o uter one. All the model s are at the same evolutionary stage, when the central H mass fraction 
is about 0.40. Figures taken from lEkstrom et~afl (|2008aT ). 
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FIG. 24 Left panels: Evolution of Z — models (with rotation: continuous lines; without rotation: dotted lines) in the 
Hertzsprung Russell diagram. The grey area shows the zone of the diagram where He burns in the core of the rotating models. 
Right panel: Evolution of the Q/Q, CI i t ratio during the MS. All the models start the M S with u eq = 800 km s^ 1 , except the 9 
M Q which starts the MS with 500 km s" 1 . Figures taken from lEkstrom et al.l (|2008aT ). 
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FIG. 25 Left panel: evolution in the Hertzsprung-Russell diagram; Right panel: evolution of the mass of the model. The 
mass indicated is the mass lost at each stage, not a summation. Black line: rotating model; continuous part: beginning of MS 
(X c = 0.753 down to 0.58; dashed part: rest of the MS; dotted part: beginning of core He-burning phase (Y c = 1.00 down to 
0.96); c ontinuous part: rest of the He-burning. Grey line: non-rotating model for comparison. Figures taken from Ekstrom ct al. 
(|2008al 1. 




B 10 20 40 60 80 100 



FIG. 26 Ranges of masses of different types of supernovae at different metallicities. The shading on the right indicates the 
area where formatio n of a black hole (BH) is expected ; elsewhere, the remnant is a neutron star (NS). Figure taken from 
iGeorgy et af] (|2009l ). 
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FIG. 27 Observed rate of type Ic SNe (triangles with error bars) and predicted rate of type Ic SNe whose progenitor is a WO 
star (dashe d line) with respect to the total number of core collapse SNe (CCSNe). The grey rectangle repre sents the extension i n 
metallicity i|Modiaz et al.1 . 120081 ) and rate l|Podsiadlowski et all 120041 ) of GRB events. Figure taken from lGeorgy et all (|2009f> . 





FIG. 28 Young bi-modal shaped nebula blown from the highly aspherical, polar dominated winds around a 20 Mq star at 
Z — 10 _J . The star has evolved during 21552 yr since it has left the ZAM S. The maximal ext ension of the wind is ~ 1.5T0 18 cm. 
This model was computed by the A-Maze code (Waldcr fc Folinil . l200Cf ). Figure taken from lGeorgy et ail (|2009l ). 
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FIG. 29 Stellar yields divided by the initial mass, pl^, as a function of the initial mass for the non-rotating (left) and rotating 
(right) models at solar metallicity. The different total yields (divided by m) are shown as piled up on top of each other and 
are not overlapping. 4 He yields are delimited by the top solid and long dashed lines (top shaded area), 12 C yields by the long 
dashed and short-long dashed lines, 16 O yields by the short-long dashed and dotted-dashed lines and the rest of metals by 
the dotted-dashed and bottom solid lines. The bottom solid line also represents the mass of the remnant {Ml^/m). The 
corresponding SN explosion type is also given. The wind contributions are superimposed on these total yields for the same 
elements between their bottom limit and the dotted line above it. Dotted areas help quantify the fraction of the total yields 
due to winds . Note t hat for 4 He, the total yields is smaller than the wind yields due to negative SN yields. Figure taken from 
iHirschi et all (|2005bh . 




FIG. 30 Product of the stellar yields, mp*^ by Salpeter's IMF (multiplied by a arbitrary constant: 1000 x Af -2 ' 35 ), as a 
function of the initial mass for the non-rotating (left) and rotating (right) models at solar metallicity. The different shaded 
areas correspond from top to bottom to mp" x 1000 x M -2 ' 35 for 4 He, 12 C, 16 Q and the res t of the heavy elements. The 
dotted areas show for 4 He, 12 C and 16 O the wind contribution. Figure taken from lHirschi et all (|2005bh . 
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FIG. 31 Yields comparison between the non - rotat ing Z = models from IChieffi fc Limongil l|2004 CL04, open pentagons), 
the rotating Z — 10~ 8 model s from IHi rschi (2007|, filled squares) and the rotating (filled triangles) and non- rotating (open 
triangle s) Z = models from lEkstrom et a l. (2008a). Left: 12 C; centre: 14 N; right: 16 0. Figure taken from lEkstrom et all 
(|2008al ). 
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FIG. 3 2 Evolution of th e N/O and C/O ratios. Data points for halo stars are from ITsraelian et al] (12004 open squares) 
and of ISpite et al] (|2005l . stars). The points with error bars are for DLA systems from iPettini et al] l|2008] ). The (black) 
continuous curve is the c hemi cal evolution model o btained with the stellar yields of slow rotating Z — 10~ 5 models from 
iMevnet fc Maed"erl ll2002af) and lHirschi et all (|2005bf ). The (magenta) dashed line includes the yields of fast rotating Z — 10~ 8 
models from IHir schi ( 200_3) at very low metallicity. The (blue) dotted curve is obtained using the yields of the Z = models 
presented in lEkstrom et all ([2008a) up to Z = 10~ 10 . The chemical evolution models are from IChiappini et~ai1 (|2006l ) . 
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FIG. 33 s-process distributions between 57 Fe and 138 Ba normalized to solar for the 25 M and [Fe/H]=-4 at the end of the 
convective C-burning shell. The horizontal line corresponds to the ls O overabundance in the C shell (thick line), multiplied 
and divided by two (thin lines). Isotopes of the same element are connected by a line. The cases presented are the following: 
i) case without additional 22 Ne (black triangles) and X( 22 Ne)i n i = 5.21 x 10 -5 ; it) case with additional 22 Ne (open squares, 
full squares for the s-only isotopes) and X( 22 Ne)i n i = 5.0 x 10~ 3 ; in) case with ad ditional 22 Ne (open circles, full circles for the 
s-only isotopes) and X( 22 Ne) in i = 1.0 x 10" 2 . Figure from iPignatari et al.1 l|200Sft . 
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TABLE I 


Grids of rotating 


massive 


star models (starting from 2000) 


Masses 


Z 


Initial 


Magnetic 


Ref. 


M 




rotation 


field 








km/s 


(interior) 




9- 200 





0-800 


no 


Ekstrom et al. (2008a) 


3-60 





39-1423 


no 


Ekstrom ct al. (2008b) 


9-85 


0.00000001 


800 


no 


Hirschi (2007) 


2-60 


0.00001 


0-400 


no 


Mevnet fe Maeder (2002a) 


20-60 


0.00001 


230-605 


yes 


Yoon fe Langer (2005) 


12-60 


0.00001 


- 936 


yes 


Yoon et al. (2006) 


3-60 


0.0001 


39-1017 


no 


Ekstrom ct al. (2008b) 


20-200 


0.0005 


0-800 


no 


Dccrcssin ct al. (2007b) 


20-60 


0.001 


230-605 


yes 


Yoon & Langer (2005) 


12-60 


0.001 


0-747 


yes 


Yoon ct al. (2006) 


3-60 


0.002 


32-879 


no 


Ekstrom et al. (2008b) 


12-60 


0.002 


0-653 


yes 


Yoon ct al. (2006) 


5-60 


0.0021 


0-600 


yes 


Brott ct al. (2011) 


9-60 


0.004 


0-300 


no 


Maeder fe Mevnet (2001b) 


30-120 


0.004 


300 


no 


Mevnet fe Maeder (2005) 


12-60 


0.004 


0-507 


yes 


Yoon ct al. (2006) 


5-60 


0.0047 


0-600 


yes 


Brott et al. (2011) 
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30-120 


0.008 


300 


no 


Mcynct & Macdcr (2005) 


5-60 


0.0088 


0-600 


yes 


Brott ct al. (2011) 


8-25 


0.020 


0-474 


no 


Heger & Langer (2000) 


8-25 


0.020 


200 


no 


Heeer et al. (2000) 


9-120 


0.020 


0-300 


no 


Mevnet & Maeder (2000) 


9-120 


0.020 


0-300-500 


no 


Mevnet & Maeder (2003) 


12-60 


0.020 


0-300 


no 


Hirschi et al. (2004) 


12-35 


0.020 


200 


yes & no 


Heeer et al. (2005) 


16-40 


0.020 


210-556 


yes 


Yoon et al. (2006) 


3-60 


0.020 


28-732 


no 


Ekstrom et al. (2008b) 


20,25,40,60,85,120 


0.040 


0-300 


no 


Mevnet & Maeder (2005) 



TABLE II Values of the average and largest [N/H] excesses observed for different types of stars in the Galaxy, LMC and SMC. 
The number of stars used to obtain the quoted values are indicated as well the sources for the observations. 

Types of stars [N/H] in Galaxy [N/H] in LMC [N/H] in SMC 

mean max. N* Rcf. mean max. N* Ref. mean max. N* Ref. 

O stars 0.3-0.5 0.4-0.6 3 [1] 1.0-1.2 1.3-1.5 14 [8] 

B dwarfs M < 20 M Q 0.2-0.4 0.5-0.7 22 [2,3,4,5] 0.3-0.5 0.9-1.1 35 [5] 0.5-0.8 1.3-1.5 34 [5] 

B giants, super. M < 20 M 0.4-0.6 0.7-0.9 10 [5] 0.6-0.8 1.1-1.3 21 [5] 

B giants, super. M > 20 M 0.4-0.6 0.8-1.0 38 [6,7] 

1] Villamariz fe Herrero (2005 ): Villamariz et al. (2002 ) 
2 Uics & Lambert (1992) ^ 

3i kiiian iiyyiii 

41 Morel'et al. (2008) 
51 Hunter et al.T2W9) 

6,71 (Jrowther et al. (2006) ; Searle et al. (2008) 
81 Heap et al. ILMU^) 



